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1. Introduction 

In this paper we would like to review one of outstanding open problems of quan- 
tum physics - rigorous understanding of the energy momentum spcctnim of homo- 
geneous Bose gas at zero temperature. We describe various rigorous and heuristic 
arguments about its shape. In particular, we discuss a number of versions of the so- 
called Bogoliubov approach. We use the main idea of this approach to give rigorous 
upper bounds on the energy- momentum spectrum of the Bose gas. 

There exists little rigorous work on this subject. We think that mathematicians 
avoid this topic not only because of its difficulty. Unfortunately, it is not easy to 
formulate questions in this domain that are, on one hand, physically relevant, and 
on the other hand, mathematically clean and precise. We try to ask a number 
of such questions, some of them rather ambitious, but some, perhaps, within the 
reach of present methods. We think that rigorous methods of spectral analysis and 
operator theory could be very helpful in clarifying this subject. 

1.1. Bose gas in canonical approach. One can distinguish two possible ap- 
proaches to the Bose gas at positive density: "canonical" fixing the density p - 
and "grand-canonical" - fixing the chemical potential /z. In most of our paper we 
will concentrate on the latter setting. Nevertheless, in the introduction, as well as 
in Section 2 we will stick to the canonical approach. 

We suppose that the 2-body potential of an interacting Bose gas is described by 
a real function v defined on M**, satisfying t;(x) = w(— x). We assume that ■u(x) 
decays at infinity sufficiently fast. 

A typical assumption on the potentials that we have in mind in our paper is 

w(k) > 0, ke (1.1) 

where the Fourier transform of v is given by 

i)(k) := ( t;(x)e-*''dx. (1.2) 
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Potentials satisfying (1.1) will be called repulsive. Note, however, that a large part 
of our paper does not use directly any specific assumption on the potentials. 

n-body Schrodingcr Hamiltonian of the homogeneous Bosc gas acts on the Hilbert 
space Lg((M'')") (symmetric square integrable functions on (M'')") and is described 
by the Hamiltonian 

H- = -j2l\+ Yl ^i^^-^j) (1-3) 

i=l ^"^iKj^n 

and the momentum operator 

n 

P" :=^-iVx.. 
1=1 

(H, P) is a collection of 1 + commuting sclf-adjoint operators, hence wc can 
ask about the properties of their joint spectrum, called the energy-momentum 
spectrum. 

(1.3) describes however only a finite number of particles in an infinite space. We 
would like to investigate homogeneous Bose gas at positive density. It is a little 
problematic how to model such a system. A natural solution would be restricting 
(1.3) to e.g. A = [—L/2, L/2]'^, the rf-dimcnsional cubic box of side length L, with 
Dirichlet boundary conditions. This will, however, destroy its translational invari- 
ance. Therefore, following the accepted although somewhat unphysical tradition, 
we consider the Bose gas on a torus. This means in particular that the potential v 
is replaced by 

^^W = ^ E ^'^■'^(k), (1-4) 

where k e ^Z'' is the discrete momentum variable and V = L'^ is the volume of 
the box. Note that is periodic with respect to the domain A, and f^(x) v{x) 
as L — » 00. The system on a torus is described by the Hamiltonian 

H^'"" = -Y^^Af + J2 ^"-(^i-^j) (1-5) 

i—l l<i<j<n 

acting on the space Lg(A") (symmetric square integrable functions on A"). The 
Laplacian is assumed to have periodic boundary conditions. 
Let us denote by i?^'" the ground state energy in the box: 

£;^'" :=infspi?^'", 

where spi^ denotes the spectrum of an operator K. 

The total momentum is given by the vector of operators 

n 

pL,n 

i=l 

Its joint spectrum equals ^Z'^. 

Clearly, H^''^ and P^'" commute with each other. Therefore we can define their 
joint spectrum 

O-TT 

sp(iI^'",P^'") cMx — Z'', 
Ij 

which will be called the energy-momentum spectrum, in the box. By the excitation 
spectrum in the box we will mean sp(ff^'" — £'^'", P^'"). 

For k s ^7/, we can define the Hamiltonian ff^'"'(k) to be the restriction of 
if^'" to the supspace of P^-" = k and infimum of the excitation spectrum (lES) in 
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the box as 



e-^'"(k) := inf sp(iJ^'"(k) - E"^'"') (1.6) 
By the infimum of the energy-momentum spectrum in the box we will mean E^''" + 



e 



(k). 



It is believed that the properties of the Bose gas simphfy in the thermodynamic 
hmit. It means that one should fix p > 0, take the number of particles equal to 
n = pV ^ and then pass to the limit L ^ x. Unfortunately, as far as we know, the 
Hamiltonians iJ^'" — iJ^'" do not have a limit as self-adjoint operators. One can 
hope, however, that the lES has some kind of a limit. 

Mathematically it is not obvious how to define this limit, since for finite L the 
lES is defined on the lattice^ ^Z^' and in thermodynamic limit it should be defined 
on R''. Below wo propose one of possible definitions of the lES in thermodynamic 
limit. 

For k e M'' and p > 0, we take S > and set 

e^'iKS) (1.7) 

liminf ( inf -i e^'"(k') : k' e |k - k'| < ^, p = 



^inf |e^'"(k') 



This gives a lower bound on the lES for the momenta k' in the window in the 

momentum space around k of diameter 2S. The quantity e''(k, (5) increases as S 
becomes smaller. The lES in the thermodynamic limit is defined as its supremum 
(or, equivalently, its limit) as 5 \ 0: 

e^(k) :=supe^(k,(5). (1.8) 

Under Assumption (1.1) it is easy to prove that E^'""- is finite and e^'"(0) = 
e''(0) = (see Theorem 3.1 and Proposition 3.3). 

Conjecture 1.1. We expect that for a large class of repulsive potentials the fol- 
lowing statem,ents hold true: 

(1) The function M'' 9 k i— > e''(k) e M+ is continuous. 

(2) Let k e K'^. If L ^ oo, UL ^ oo, ^ ^ p, \ll e ^'L'^, and k^ k , we 
have that e^'"^(kL) ^ e''(k). 

(3) Ifd>2, then inf ^ =: c„ > 0. 

(4) There exists the limit lirti =: Cph > 0. 

(5) The function M** 9 k i-^ e''(k) is subadditive, that is, e''(ki -|-k2) < e''(ki)-|- 

e''(k2). 

Statements (1) and (2) can be interpreted as some kind of a "spectral thermo- 
dynamic limit in the canonical approach". Note that if (1) and (2) are true around 
k = 0, then we can say that there is "no gap in the excitation spectrum" . 

The properties (3) and (4) of the Bose gas were predicted by Landau in the 
40's. Shortly thereafter, they were derived by a somewhat heuristic argument by 
Bogoliubov [3]. 

(3) is commonly believed to be responsible for the superfluidity of the Bose gas. 
More precisely, it is argued that because of (3) a drop of Bose gas travelling at 
speed less than Ccr will experience no friction. This argument is described e.g. in 
the course by Landau-Lifshitz [14] and in [39], see also Section 2.6. 

Note that in dimension d = 1 the statement (3) should be replaced by 

(3)' If d= 1, then 

e''(k-F27rp) =e^(k). (1.9) 



The statement (3)' has a simple rigorous proof, which we wiU describe in our 
paper. It impUes that in dimension d = I the excitation spectrum is periodic with 
the period 2T:p. It follows by the well known argument that involves boosting all 
particles simultaneously by the velocity ^ (see e.g. [17]). 

Excitation spectrum with the property described by (4) is often described as 
phononic and the excitations with such a spectrum are called phonons. One also 
expects that the parameter Cph coincides with the speed of sound - a parameter in 
principle macroscopically measurable. 

Let us describe a heuristic argument for (5). Suppose that excitations can be 
described by certain elementary quasiparticles with a dispersion relation k i— > ix'(k). 
We assume that the state consisting of quasiparticles with momenta ki , . . . , k„ has 
the excitation energy a;(ki) + • • • + tj(k„). Then it is easy to sec that the lES is 
the subadditive hull of u!(k), that is 

e(k) =inf{cc;(ki) + ---+w(k„) : ki + • • • + k„ = k, n = l,2,...}, (1.10) 

which is the largest subadditive function less than e. In Appendix B we describe a 
somewhat more elaborate, but still heuristic, argument that seems to indicate that, 
in thermodynamic limit, the Bose gas has a subadditive excitation spectrum. 

Note that free Bose gas does not satisfy Conjecture 1.1. In this case a;(k) = k^ 
and e(k) = 0, see Fig. 1.1 
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Figure 1.1. Excitation spectrum of the free Bose gas 

There are not so many subadditive functions. There exist, however, subadditive 
functions, which satisfy the properties described in our conjecture in (3) or (3)', 
and in (4). We discuss basic properties of subadditive functions in Appendix A. 
We could not find these facts in the literature, although they probably belong to 
the folk knowledge. 

1.2. Experimental evidence. To our experience, most physicists interested in 
this subject (but not all) would agree that one should expect Conjecture 1.1 (as 
well as the analogous Conjecture 3.4 formulated in the grand-canonical setting) 
to be true. Let us start with a brief account of experimental evidence for these 
conjectures. 

Theoretically, the cross-section for neutron scattering against a droplet of Helium 
IV at zero temperature is approximately proportional to the so-called van Hove 
formfactor S'(w,k) ([36], see also (H.6)). S'(ti;,k) is a measure of the density of 
excitations of the Bose gas at energy lo and momentum k at zero temperature. 
Therefore, ^(w, k) is zero below the curve k i— > e''(k). It seems reasonable to 
suppose that more is true: S{uj,k) should be nonzero everywhere above the curve 
k H- » e''(k). If in addition Conjecture 1.1 is true, then the lower boundary of the 
support of S{uj,'k) should satisfy (3), (4) and (5) of this conjecture. 
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To our understanding, within experimental accuracy, experiments on Helium IV 
at low temperatures seem to confirm the above theoretical expectations. 

Actually, experiments seem to say more than this. At least for low momenta, one 
observes a sharp peak of 5(0-1, k) along a curve similar to k 1-^ "-^(k) at Fig. 1.2, see 
[38] and Fig 1 of [22] . This curve is interpreted as the dispersion relation of a quasi- 
particle (elementary excitation spectrum). These quasiparticles are called phonons 
for small momenta and rotons around the local minimum of the dispersion relation. 
To our understanding, experiments indicate that below the subadditive hull of the 
elementary excitation spectrum the value of S{u^ k) drops down substantially. (In 
the case of Fig. 1.2, this subadditive hull equals k 1-^ e(k)). 

Experiments involving excited phonon states are usually successfully interpreted 
in terms of multi-quasiparticle states whose momenta and energies are additive [22] . 
This also implies that the energy of multi-quasiparticle states lies above k 1-^ e(k). 

A similar picture arises in the case of Bose-Einstein (BE) condensates of al- 
calic metals. For example, the reader can consult Fig. 2 of [31], which shows the 
quasiparticle spectrum of the BE condensate of *^Rb around zero temperature. 
Compared to Helium IV, the main difference is the absence of the rotonic part of 
the elementary excitation spectrum, see Fig. 1.3. 

Of course, it is difficult to interpret real experiments in terms of rigorous state- 
ments. The setup that we describe in this paper does not apply in all its details 
to realistic BE condensates. First of all, both in Helium IV and alcalic metals, 
the potential has typically an attractive part and a hard core. Therefore, strictly 
speaking it does not belong to the class that we would like to consider in this paper. 
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In the case of Helium IV, the situation is further complicated by the fact that 
the Schrodinger operator of the form (1.3) is not believed to describe it adequately 

- 3-body interactions are probably relevant. This problem does not appear in BE 
condensates of alcalic metals, where apparently one can assume that only 2-body 
interactions play the role. 

BE condensates of alcalic metals have a different conceptual problem absent in 
the case of Helium IV: they do not represent the true ground state but only a 
metastable state. 

1.3. Bogoliubov approximation. There are many theoretical physics papers de- 
voted to the Bosc gas. To our surprise, their authors usually avoid making precise 
statements or conjectures about the excitation spectrum of the Bose gas. (A notable 
exception is [39], where a definition of the lES similar to (1.7) can be found). 

In [3] Bogoliubov proposed an approximation, which implies that the Bose gas 
should be described by elementary excitations with the spectrum 



Within this approximation, the lES equals e(^g, the subadditive hull of (see 
(1.10)), which has the properties described in Conjecture 1.1. 

Note that if we replace the potential v with Xv and the density p with p/X (with a 
positive A), then neither u>^^ nor e^^ depend on A. In fact, it is natural to conjecture 
that e(Jg describes the true lES in the weak coupling/large density limit. 

More precisely, let e'''''*(k) be the lES for the potential Aw. 

Conjecture 1.2. Let d>2. Then for a large class of repulsive potentials we have 



Note that Conjecture 1.2 is certainly wrong in dimension d = 1, because of (1.9). 

We do not know complete proofs of a statement similar to Conjectures 1.1, 1.2, 
as well as their grand-canonical analogs described later on in our paper. We believe 
that to prove or disprove them would be an interesting subject for research in 
mathematical physics. 

Many theoretical works on the energy-momentum spectrum of the Bose gas 
instead of the correct Hamiltonian iJ^-" (or its second-quantized version and 
the grand-canonical version -ff^) consider its modifications. They either replace 
the zero mode by a c-numbcr or drop some of the terms, or do both modifications 
[1, 13, 8, 10, 34, 40]. These Hamiltonians have no independent justification apart 
from being approximations to the correct Hamiltonian in some uncontrolled way. 
Let us stress that in our paper we are mostly interested in the correct Hamiltonian 
and not its modifications: all our statements will be related cither to or to 

the grand-canonical Hamiltonian (the second quantization of H^ " — fm). 

1.4. Organization of the paper. The paper is divided into several sections and 
appendices. The individual sections use sometimes slightly different notation and 
arc devoted to different apsects of the Bosc gas. 

Section 2 is devoted to some facts about Bose gas that are naturally formulated 
in the canonical approach, where we start with a definite number of particles and 
go to thermodynamic limit keeping the density fixed. The later part of the paper, 
where we use the grand-canonical approach, is independent of this section. 

We start with a discussion of the Galileian covariance in a finite box with periodic 
boundary consitions. We believe that this is relevant if one wants to understand 




(1.11) 



lime^'^(k)=e^g(k). 
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the physics of the Bose gas. Even though what we present is elementary, we did 
not find most of it in the literature. 

In Subsection 2.4 we discuss the case of dimension d = 1. We prove that the 
excitation spectrum is periodic in momentum. This fact is known to some experts 
and we do not claim its discovery - nevertheless, we have never seen it explicitly 
stated in the literature. 

In Subsection 2.6 we describe an argument that links the properties of the exci- 
tation spectrum to superfluid behavior. The argument that we describe is slightly 
different from the one usually stated in the literature [18, 39] - in particular it 
applies to systems confined to a finite volume. 

One of the most imporant quantities in in superfluidity is the so-called critical 
velocity. There are several non-equivalent definitions of this concept in thermody- 
namic limit. Wc argue, that to obtain the definition that is relevant for superfluidity 
one should first take thermodynamic limit for the excitation spectrum, and only 
then compute the critical velocity. 

In Subsection 2.7 wc present the variational ansatz due to Bijls [2] and Feynmann 
[7] for the excitation spectrum. 

In Subsection 2.9 we describe the (non-rigorous but interesting) argument due 
to Onsager [25] that indicates the phononic character of the excitation spectrum 
obtained by this ansatz. Our presentation follows that of a recent paper [39]. 

Section 3 is the central part of our paper. Starting with this section, we switch 
to the grand-canonical setting. This means that we allow the number of particles 
to vary and we fix the chemical potential fi. We also use the formalism of second 
quantization. 

In Subsection 3.1 we describe the formalism and formulate Conjecture 3.4, the 
grand-canonical analog of Conjecture 1.1. 

In Subsection 3.2, we discuss the Hamiltonian obtained by a c- number substitu- 
tion of the zero mode. We describe the theorem of Lieb, Seiringer and Yngvason 
[20] saying that this approximate Hamiltonian gives the correct energy density in 
thermodynamic limit. Note that the result of [20] is more general, it concerns an ar- 
bitrary temperature. Our presentation sticks to the zero temperature, which allows 
for some minor simplifications. 

In Subsection 3.3 we describe the Bogoliubov approximation. Its original form 
was formulated in the canonical setting of fixed density in the second-quantized for- 
malism. We follow its grand-canonical version, which can be traced back to Beliaev 
[1] and Hugenholz - Pines [13], see also a recent review paper by Zagrebnov and Bru 
[40] . Even though this is a classic reasoning, our presentation seems to be somewhat 
different from and more satisfactory than what we have seen in the literature. Its 
first step is a variational problem involving coherent states. The second step is the 
Bogoliubov translation and rotation adapted to the resulting approximate ground 
state. Our reasoning does not involve the c-number substitution: we treat the zero 
mode quantum mechanically. 

One can try to improve on Bogoliubov's approximation by looking for the mini- 
mum of the energy among translation invariant squeezed states. To our knowledge, 
in the context of the Bose gas this idea first appeared in the paper of Robinson 
[30]. Robinson considered a slightly more general class of states - quasi- free states. 
He noticed, however, that in the case he looked at it is sufficient to restrict to pure 
quasi-free states, which coincide with squeezed states. One should mention also [6], 
where a variational bound on the pressure of Bose gas in a positive temperature is 
derived by using quasi-free states. In the literature this approach often goes under 
the name of the Hartree-Fock-Bogoliubov method. 
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Only an upper bound to the ground state energy is considered in [30] . We go one 
step further: we show how this method can be extended to obtain upper bounds on 
the infimum of the energy-momentum spectrum by using one-particle excitations 
over squeezed states. 

After finding the minimizing squeezed vector ^, it is natural to express the 
Hamiltonian in the new creation/annihilation operators b^^/bk, for which the new 
approximate ground state is a Fock vacuum. We show that the resulting quadratic 
Hamiltonian has no terms involving 6, b*, bb and b*b*. The Hamiltonian becomes 

H = C + ^ £'(k)&k^k + terms of order 3 and 4. (1.12) 

k 

This new form of the Hamiltonian yields immediately an interesting estimate 
on the energy-momentum spectrum. In fact, the vectors have precisely the 
momentum k and the energy equal to C + D{k). One could ask how good is this 
estimate. It turns out that it has a serious drawback. As we show in Subsection 
3.7, under quite general circumstances wo have D{0) > 0. 

Perhaps, this is the most important (even if negative) finding of our paper. It 
implies that at the bottom of its spectrum 1-particle excitations over a squeezed 
state arc poor test functions for the excitation spectrum. 

In the literature, the existence^ of a gap in various approximation schemes that 
try to improve on the original Bogoliubov's one has been noticed by a number of 
authors [10, 34]. However, to our knowledge those authors did not consider the 
correct Hamiltonian, but always used one of its distorted versions. 

The Hilbert space of the homogeneous Bose gas can be naturally factorized into 
the tc;nsor product of an infinite family of Hilbert spaces for various values of the 
momenta. Variational ansatzes involving translation invariant squeezed states, as 
well as particle excitations over the squeezed states have a common feature - they 
are factorized with respect to this tensor product. We call such states uncorrelated. 
One can pose a question: how good are uncorrelated states as variational test 
functions in many body problems? It seems to us that they have serious drawbacks 
- in particular, we conjecture that they typically give spectrum with an energy gap. 

In Section 4 we discuss approaches to the Bose gas based on perturbation theory. 
We would like to treat the coupling constant A as a small parameter, keeping the 
chemical potential /i fixed. 

In Subsection 4.1 we describe a naive splitting of the Hamiltonian into a main 
part and a perturbation based on the usual Bogoliubov approach. Unfortunately, 
this approach seems to fail because of a serious infrared problem. We also formu- 
late Conjecture 4.1, which is the grand-canonical analog of Conjecture 1.2 and is 
suggested by Bogoliubov's approximation. 

In Subsection 4.2 we propose a certain systematic procedure for perturbation 
expansion, which avoids the infrared problem. This procedure uses (1.12) as the 
starting point for the expansion. The 3rd and 4th order term are treated as per- 
turbations. The advantage of this procedure is that it does not drop any terms 
from the Hamiltonian. All the works on the Bose gas based on perturbation the- 
ory [1, 13, 8] that we know involve the c-number substitution. This substitution, 
even if justified for the energy density [20], is unfounded for finer quantities such 
as the infimum of the excitation spectrum. Our perturbative procedure does not 
involve distorting the Hamiltonian. Therefore, in our opinion, it is superior from 
the physical point of view. 

In Section 5 we describe various inequalities on the Bose gas that can be proved 
rigorously. These results are consistent with the absence of the energy gap and the 
phononic form of the excitation spectrum. They are obtained by relatively simple 
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methods, involving especially the so-called /-sum rule. Our presentation is based 
on the work of Bogoliubov [4] and on results presented by Stringari [32, 33]. 

In appendices we present some background material, to make our paper accessi- 
ble to a larger audience. In Appendix F we describe technical computations. 

1.5. Additional remarks about the literature. Let us make some additional 
remarks about the literature of the subject. The case of dimension d = 1 and 
repulsive delta interactions has been studied in detail. Girardeau [9] studied the 
case of "infinite" coupling (which amounts to the Dirichlct boundary conditions). 
The case of an arbitrary positive coupling constant was studied in [18, 17], where 
arguments for the absence of a gap and the phononic shape of the excitation spec- 
trum in the thermodynamic limit arc given. [39] gives a full rigorous proof for the 
linearity of the excitation spectrum in thermodynamic limit for Girardeau's model. 
Note, however, that the 1-dimensional case is believed to be quite different from 
the case d> 2. 

There exists a large literature on the energy density of the Bose gas. The energy 
density can be defined as 

pf := lim -, (1.13) 

where p = y is kept fixed. There exist derivations of the asymptotics of e'' in 
dimension d = 3 for small p going back to [5, 13, 15, 16]. This asymptotics is not 
restricted to small potentials - it covers also the case of hard-core potentials. One 
can show rigorously (see [21] and references therein) that the leading term of this 
asymptotics correctly describes the energy density. Similar results can be shown in 
dimension d = 2. 

Note that the energy density is easier to study than the infimum of the excitation 
spectrum. Besides, it does not capture some interesting physical phenomena the 
excitation spectrum is responsible for. Note also that the above mentioned results 
involve the following limit: the scattering length is kept fixed (this can be achieved 
e.g. by fixing the potential) and the density goes to zero. In our paper we usually 
consider a different limit: the chemical potential is kept fixed and the coupling 
constant in front of the potential goes to zero. One can have various opinions 
comparing the physical relevance of the two limits. In any case, the latter limit 
seems more appropriate if one wants to capture the phononic character of the 
excitation spectrum. 

Another direction of rigorous research involves studying the so-called Gross- 
Pitaevski limit. Again it concerns mostly the dimension d = 3. The quantity 
that is kept fixed is where a is the scattering length and n goes to infinity. 
The Gross-Pitaevski limit is usually presented with a fixed L and the scattering 
length a going to zero, which is achieved by an appropriate scaling of the potential. 
Equivalently, one can fix the potential, consider L ^ oo. and scale the density as 
p ~ L^""^ as L — > 00. In this limit, Lieb, Seiringer and Yngvason have obtained 
a number of interesting and precise results [21]. In particular, they are able to 
approximate the behavior of the n-body system by a non-linear effective equation 
- the Gross-Pitaevski equation. Note, however, that in this limit it is difficult to 
say something interesting about the excitation spectrum, because the density, and 
hence the speed of sound and the critical velocity, go to zero. 

The result of [19] (see also Theorem 5.3 of [21]) can be interpreted as the posi- 
tivity of the critical velocity in finite volmnc. The critical velocity, as understood 
in this result, goes to zero in thermodynamic limit. It is presumably far too low to 
fully account for the physical phenomenon of superfluidity, see Subsection 2.5 for a 
discussion. 
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Finally, let us note that there exists a number of interesting rigorous results 
about Bose gas in positive temperatures [4, 6, 12, 32, 33, 26, 20]. 
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2. Canonical approach 

In this section we will always work on a Hilbert space of fixed number of particles. 
We will use the Hamiltonian iJ^'" introduced in the introduction. To simplify the 
notation, we usually drop the superscripts n,L, writing e.g. H instead of H^'". 

2.1. Free Bose gas. In finite volume, the momentum is restricted to 

k = —k, k e Z"^. 

L 

It is easy to compute exactly the excitation spectrum of the free Bose gas in a finite 
volume, see Fig 2.1. In particular, in dimension c? = 1, its infimum is given by the 
broken line with vertices at 

k^^k, e^.«(k) = ik2, kez. 

In an arbitrary dimension, we just add the contributions from each dimension: 
e(ki,...,k<j) ==e(ki) + --- + e(kd). 
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Figure 2.1. Infimum of excitation spectrum of free Bose gas in 
finite volume 



2.2. Galileian covariance on a box vifith periodic boundary conditions. 

In infinite volume the Galileian covariance involves an arbitrary value of velocity. 
This is not the case on a box with periodic boundary conditions (torus), where the 
Galileian covariance is only rudimentary. To describe it we will restrict ourselves to 
boosts in the first coordinate. The following operator, which we will call the boost 
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operator adds simultaneously velocity ^ to all particles in the direction of the first 
coordinate: 

Ui exp I — 2_,^»i I • 

(xii denotes the 1st coordinate of the ith particle). Clearly, Ui preserves the domain 
of H and is a unitary operator on Ls{A) satisfying 

2nn 
IT' 



UiPiU* 



Pi 



2n (27r)2n 

(Pi denotes the first component of the total momentum). Hence 



(2.1) 
(2.2) 



UiiH- —p^ \u: = H - —P^ 

^ 2n J ^ 2n 



(2.3) 



(2.3) and (2.2) imposes a severe restriction on the shape of the excitation spectrum: 
has to be invariant with respect to translations by see Fig. 2.2. 
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Figure 2.2. Typical infimum of excitation spectrum of interacting 
Bose gas in finite volume 



2.3. Critical velocity. The critical velocity in finite volume is defined by 

. 6^'"(k) 

C(,; := mf — — — . 

k#0 |k| 

Recall that we are interested in thermodynamic limit, which involves L, n 



with 



p>0. 



Note that for the free Bose gas we have 



L.n 



Thus the critical velocity is positive, but goes to zero in thermodynamic limit. 

In the interacting case we have a sequence of low energy states with the momen- 
tum and the excitation energy obtained by boosting the ground state: 

2^~ ,_(M!„k2 



where k £ Z"*. Expressed in terms of density this gives 

(2^)^ , 



k = 27r/9L'^"^k, e 



Therefore, in the general case the critical velocity is not greater than in the case of 
the free gas: 



< 



In dimension > 2, the momentum of these states escapes to infinity, so in a 
sense they are not visible in thermodynamic limit. 

2.4. Bose gas in dimension d — 1. Bosonic gas in dimension d = 1 seems to 
have different properties than in higher dimensions. In particular, statement (3) of 
Conjecture 1.1 should be replaced by (3)', that is 

Theorem 2.1. M 9 k i-^ ^''(k) is periodic with the period 27rp. 




Figure 2.3. Typical infimum of excitation spectrum of Bose gas 
in dimension d = I 

Theorem 2.1 easily follows from the invariance of the spectrum (2.4). 

2.5. Twisted boundary conditions. One can also consider the boost operator 
for an arbitrary velocity. In fact, for a G R define the unitary operator 



Ui{a) := cxp [ 



Let [a] denote a mod {2tt). We will write -PiJq], Lf[a] for the momentum and the 
Hamiltonian with the boundary condition in the 1st coordinate twisted by [a] . (It 
means that the elements of the domain of these operators for j — 1, . . . ,n satisfy 



$ ( xi, . . . ,x,- - ^Lei, 



e'" $ XI 



1 



where ei is the unit vector in the direction of the Is coordinate). 

For a not equal to a multiple of 27r the operator Ui{a) does not preserve the 
domain of H and P. Instead of identities (2.1) and (2.2) we have 

Uiia)PiUii-a) = Pi[„j-^, 



(2.5) and (2.6) imply 

= C/i(a) (h + ^Fi) U,{-a) + 



(2.5) 
(2.6) 

(2.7) 
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In particular, if < ao < tt, then 




if and only if 



H[a] - E> — , \a\ < ao. 



(2.8) 



(2.8) was used by Lieb, Seiringer and Yngvason as a criterion for the positivity of 
critical velocity in [19], see also Theorem 5.3 of [21]. 

2.6. Superfluidity. Let us describe one of experiments that show superfluid prop- 
erties of the Bose gas. 

A laser beam playing the role of an optical spoon [29] is directed into a sample of 
a Bose gas at a sufficiently low temperature. The beam makes a rotating motion. 
If the velocity of this motion is lower than a certain critical value, then the sample 
heats up very slowly. For velocities above this critical value, the sample heats up 
much faster. 

Let us try to describe an idealized mathematical model of this phenomenon, 
which is a version of the well known argument due to Landau described e.g. in 



Since our Bose gas has periodic boundary conditions, we can make an idealized 
assumption that the "laser beam" travels forever with a constant velocity w = 
{w,0, . . . , 0) G R'^. We will model it with a weak travelling potential u(x — tw) 
interacting with all particles. Thus the system is described by the Schrodinger 
equation with a time-dependent Hamiltonian 



[14, 39]. 




(2.9) 



Let us replace with 

^'t(xi, . . . ,x„) ^'t(xi 4- iw, . . . ,x„ + tw). 
We obtain a Schrodinger equation with a time-independent Hamiltonian: 





k 



Figure 2.4. Infimum of excitation spectrum of tilted Hamiltonian 
for velocity below Ccr 



14 




We know that the operator H has a nondegenerate ground state energy E. 
The corresponding ground state 4* is therefore stable with respect to smah time- 
independent perturbations. We ask the question whether the state ^E* is stable 
against a small travelling perturbation of the form (2.9). 

Let us first describe a slightly dishonest version of the argument for superfluidity 
(which actually is usually found in the literature). Let Ccr denote the critical ve- 
locity. For |w| < Ccr, the excitation spectrum of the "tilted Hamiltonian" H — wP 
looks as at Fig. 2.4, so 5* is its ground state. Hence ^E" hence is stable. For |w| > Ccr 
the excitation spectrum looks as at Fig. 2.5. Therefore ^' is not a ground state 
oi H — wP and its energy is close to energies of many other states. Hence 4* is 
unstable. 

Note that we cheated a little in the above argument. The situation that we 
described involved a finite volume, but the pictures were (as we believe) typical for 
the thermodynamic limit. The actual plot in finite volume resembles Figure 2.6. 
In particular the critical velocity in finite volume c";^ is small and goes to zero 



-V-l- 



4H- 



2jw 



2ir« 



Figure 2.6. Infimum of excitation spectrum of "tilted Hamilton- 
ian" in finite volume 

15 



in thermodynamic limit. Physical evidence seems to indicate that the observed 
phenomenon of superfluidity involves much larger velocities, which are positive 
in thermodynamic limit. Note that the Figs 2.4 and 2.5 represent the excitation 
spectrum in thermodynamic limit and do not show the low lying states present in 
flnite volume, which have a very low critical velocity. 

Let us try to present a more physical argument for superfluidity. Suppose that 
the Fourier transform of u is supported in the ball |k| < R. Define the critical 
velocity below the momentum R as 

Cer,fl:=inf|4^ k^O, |k|<i?l. (2.10) 



If the "tilted Hamiltonian" H — wP has no other eigenstates of energy close to E 
and momentum less that R, then the state ^ will be stable against the perturbation 
(2.9), at least in the 1st order. This is the case for |w| < Ccr.B- On the other hand, 
for |w| > CcT,R, we can expect many states with energy close to E and momentum 
less than R, so the stability will be lost. 

Heuristic arguments (e.g. the Bogoliubov method described later on) suggest 
that in dimension d> 2 the critical velocity c^^. is bounded away from zero even 
in thermodynamic limit and for an arbitrary R. This can be formulated in the 
following conjecture: 

Conjecture 2.2. Fix density p. Then 

p 1- L.n ^ 

' L— »oo ' V 

exists, where c^j.'*^ is defined as in (2.10). Moreover, in dimension d>2, 



inf <«>0. (2.11) 



Note in parenthesis, that Conjecture 2.2 is stronger than Conjecture 1.1 (3). In 
fact, the left hand side of (2.11) is less than or equal to c^j. introduced in Conjecture 
1.1 (3). 

2.7. Bijls— Feynmann's ansatz. For k G set 

n 

iVk-^e"'''-, (2.12) 

i=l 

acting on Lg(A"). 

It is well known that the ground state of H is nondegenerate. Denote it by 
We will write (see Appendix G) 

(A) := (4-1^*), {{A, B)) := {A{H - E)-^B) + {B{H - E)-^A). 

Note the following identity: 

i[iV^,[lf,iVk]] = yn. (2.13) 
It implies the so called f-sum rule: 

liN^iH - E)Nu) + ^(iVk(if - E)NZ) = ^n. (2.14) 
By the reality of * and H, (2.14) also equals {N^{H - E)Ni,). Define 

Sk := n-\N^Nu); (2.15) 
Xk := n-'{{N^,Ni,)). (2.16) 
By the reality of * and H, (2.16) also equals 2n-^{N^{H - Sy^Ni,). 
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By the Schwarz inequality, wc obtain 

Sk<^|k|v^. (2.17) 

Bijls [2] and Feynmann [7] proposed to consider the following variational ansatz: 

*k :-iVk*/||iVk*||, 
to obtain the excitation spectrum of the Bose gas. 
Theorem 2.3. We have 

P*k = k*k, (2.18) 

^-0 = (2.19) 

{^^\{H-E)^^) = M!>i. (2.20) 
Proof. To see (2.20) we note that 

{N^{H - E)N^) + {Nu{H - E)N^) 



{^k\{H-E)^k) = 



2(iV*iVk) 

k2 11.1 



> 



2sk A/Xk 

In the above we apply the /-sum rule to the numerator and used (2.17). □ 
It is believed that the Bijls-Feynman ansatz gives the correct behavior of the 
excitation spectrum for low momenta, and in particular it gives the value of Cph, 
which was defined in Conjecture 2.2 (4) [33]. Let us formulate this as a conjecture: 

Conjecture 2.4. Fix density p. Let R'^ 9 k i-^- egp(k) be defined as in (1.7) and 
(1.8) with (*k|(if^'" - -B^'")*k) replacing e^'"(k). (In other words, e^p(k) is a 
thermodynamic limit of the left hand side of (2.20)). Then 

e''(k) e'^p(k) 
lim —V- = lim • (2.21) 

k^O |k| k^O |k| 

2.8. Speed of sound. Recall that under broad conditions wc arc able to prove 
the existence of the energy density (1.13), which for typographical reasons in this 
subsection we will be written e(p) instead of e^. 

Another important macroscopic parameter, which in principle can be measured 
experimentally, is the speed of sound, denoted Cg. It is related to the energy density 

by 

c, = Vpe"(p) (2.22) 

(see e.g. Appendix C). 

It is believed that for low momenta 

lim #^ = T^, (2.23) 
k-^o |k| 2cs ' 

and hence {'^y\{H — i?)5'k) ^ Cs|k|. If this is the case and if the speed of soimd 
is nonzero, then the excitation spectrum given by the Bijls-Feynman ansatz is 
phononic, that is the limit on the right hand side of (2.21) exists and equals Cg. 

2.9. Relation between the speed of sound cuid Xk- Instead of arguing for 
(2.23), it seems easier to justify the relation 

limxk = l (2.24) 

In this subsection we describe a series of heuristic arguments indicating that (2.24) 
holds in thermodynamic limit. Given (2.24) and if the speed of sound is nonzero, 
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(2.20) gives a phononic lower bound on (^ftK-ff — -E)\l/k;)- This inequality is attrib- 
uted to Onsager [25]. Our presentation follows that of [39]. 

For a given vector ^ of norm one its 1-particle density is defined as 



(x) := n J ■ ■ ■ J X2, . . . x„)$(x, X2, . . . , x„)dx2, . . . dx„. 



If $ = ^f, the ground state of H, then is a constant equal to p. 
The following claim seems to be widely accepted: 

Conjecture 2.5. If L is large and ^ is a slowly varying vector close in some sense 
to the ground state then 

« J e(p$(x))dx. (2.25) 

The next claim follows from (2.25) by an application of a perturbation argument: 

Conjecture 2.6. Let u be a periodic function on R'' (or its restriction to A, where 
we assume that L is a multiple of the period of u). Assume that the mean value of 
u is zero and that u varies slowly. Then for large L we have 

^ ^u(x)2dx = ^lYu{^i){H-EY^Yu{^i)). (2.26) 



Ve"{p) J ' ' ' ' ^ 

Let us give an argument in favor of Conjecture 2.6. Consider the perturbed 
Hamiltonian 

Hr:=H + TY,<^i)- (2-27) 

i 

Let Et be the ground state energy of Ht. By the usual perturbation theory, 

^{Er-E) = 0, (2.28) 

^{Er-E) = -2l^u{^i){H-E)-'Y,<^)y (2-29) 

Applying the principle (2.25), we obtain 

Er = inf |($|ir$) + T /p<i>(x)u(x)dx : ||$|| = 1 



w i^fj/ e(pi(x))dx + T j pi(x)u(x)dx| (2.30) 

« Ve{p) 

+ inf|y' f^(pi(x)-p)2dx + r|(pi(x)-p)u(x)dx|, (2.31) 

where in (2.30) and (2.31) wc minimize over positive functions pi such that J pi (x)dx : 
n. The minimum of (2.31) is attained at 

Therefore, 

^2 



Er « l/e(p)-/^.(x)^dx. 



(P) 

Invoking (2.29) we obtain (2.26). 

Now (2.24) will follow from (2.22) and the following claim: 
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Conjecture 2.7. Xk is well defined in thermodynamic limit and satisfies 

lim Yk = ,} , ■ (2.32) 
k^o^*^ e"{p)p ^ ' 

Let us justify (2.32). We assume that L is large. Clearfy, 

+ ^^sin(kx,)(i? - i?)-i ^sin(kx,)^ 

+2Im sin(kxi)(H - E)'^ ^ cos(kxj)^ . (2.33) 

The ground state is a real vector. The Hamiltonian H, and hence also {H—E)~^, 
is a real operator. Therefore, the last term in (2.33) is zero. 

We note that coskx and sinkx are slowly varying for small k. Therefore, we 
can apply (2.26) and obtain 



(7V^(ff-E)-i7Vk) « :j;^|(cos2(kx)+sin2(kx))dx 



2 
V 

1 



e"{p) 



3. BOGOLIUBOV APPROACH IN THE GRAND-CANONICAL SETTING 

3.1. Grand-canonical approach to the Bose gas. As realized by Bogoliubov 
[3], even if one is interested in properties of the Bose gas with a fixed but large 
number of particles, it is convenient to use the second quantized description of the 
system, allowing for an arbitrary number of particles. 

An additional reformulation of the problem was noted already by Beliaev [1], 
Hugenholz - Pines [13] and others. Instead of studying the Bose gas in the canonical 
formalism, fixing the density, it is mathematically more convenient to use the grand- 
canonical formalism and fix the chemical potential. Then one can pass from the 
chemical potential to the density by the Legendre transformation. 

More precisely, for a given chemical potential > on the symmetric Fock space 

r,(L2(A)):= ® Ll{A"), 

n=0 

we define the grand-canonical Hamiltonian 

ff^ := ® (iJ^'"-Mn) 



n=0 



I 



a* (-^Ax - A«)axdx 



+ \ J J axay'y^(x- y)ayaxdxdy. 
The second quantized momentum and number operators are defined as 

T OO f 

:= n = I axflxdx, 

«=o ./ 



g pn,L ^ _j /"a*v^axdx. 

n=0 J 
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It is convenient to pass to the momentuni representation: 

k 

+ ^ S{ki+k2-k3-k4)v{k2-'k3)al^al^ak^ak^, 

ki,k2,k3,k4 

= $^<ak, 

k 
k 

where we used (1.4) to replace v^{x) with the Fourier coefRcients 'O(k). Note that 
v(k) = v{-k), and ^ T^-^/^ e"''' Ok- 

The ground state energy in the grand-canonical approach is defined as 

= inf spF^ = iiif (E^>" - /xn). (3.1) 

is a decreasing concave function of /x. To go back to the canonical condition 
(fixed number of particles) one uses 

n = -d^Ef:. (3.2) 

Both E^'"^ and Ef^ are finite for a large class of potentials, which follows from a 
simple and well-known rigorous result, which we state below. 

Theorem 3.1. Suppose that v{k) > 0, v{0) > and v{0) < 00. Then H^''" and 
are bounded from below and 

E--- > #"^-^n, (3.3) 

Proof. Let us drop the subscript L. Set 

:= 13«k+q«k= /e'^i'aiaxdx, (3.5) 
k 

(which is the second quantized version of (2.12)). Then, by a simple commutation, 
using that Nq = N and that v{k) is non- negative, we obtain 



Hij, > -jiN + ^ E 5{ki + k2 - k3 - ki)v{k2 - ks)al,^al,^a\^^a\^^ 

ki ,k2,k3,k4 

= ^E^(k)^kA^k-(M + ^)^ 



2^ k 

> + (3.6) 

Setting /i = and = n in (3.6) we obtain (3.3). Minimizing (3.6) over A'' proves 

(3.4). □ 

For k g ^Z'' let H^ik) denote the Hamiltonian Hj^ restricted to the space 
P = k. The lES in the box is efined as 

e^(k):=infspiJ^(k). (3.7) 
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For k e M'' we define the lES at the thermodynamic limit 



e„(k) := sup liminf inf efAk') . (3.8) 

Throughout most of our paper, the chemical potential n is considered to be the 
natural parameter of our problem. It often can be proven that the energy density 
exists in thermodynamic limit for a fixed /U > 

■■= lim -jf . (3.9) 

and is related to of (1.13) by the Legendre transformation 

= inf{e'' - fip}. 
p 

By definition, is decreasing and concave. Hence it is diffcrcntiablc almost every- 
where. At the points of differentiability, we can pass from the grandcanonical to 
the canonical approach by 

-d^,e^ = p (3.10) 

At the points where (3.10) has a unique solution i^{p), we should be able to relate 
the canonical and the grandcanonical lES: 

Conjecture 3.2. e''(k) = e^(p)(k). 

The following proposition is one of few rigorous facts that can be easily shown 

about the lES 

Proposition 3.3. At zero total momentum,, the excitation spectrum has a global 
minimum where it equals zero: e^'"(0) = e''(0) — and e^(0) = £^(0) — 0. 

Proof. Each iJ^^" is a non-degenerate eigenvalue of iJ^'", and commutes 
with the total momentum and space inversion. Thus each iJ^^" corresponds to zero 
total momentum, and hence by (3.1) so does E^. Hence e^'"(0) = e^(0) =0. □ 
Let us now formulate the conjectures about e^(k) (analogous to the Conjecture 
1.1 about e''(k)): 

Conjecture 3.4. We expect the following statements to hold true: 

(1) The map M'' 9 k ^/^(k) S R-i- is continuous. 

(2) Let k e K'^ If L ^ oo, kL e ^Z'^, and kz, k, we have that e^(kL) 

(3) Ifd>2, then infk^o ^ =: c„ > 0. 

(4) The limit limk^o =• ^ph > exists. 

(5) k>-^ e^{k) is subadditive. 

3.2. c-number substitution. One of the steps of the Bogoliubov method consists 
in replacing the operators Qq, ao with c- numbers: 

cto^o ~ |ck|^j 0-0 ~ a, Oq « a. (3-11) 
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This means replacing by 



2^ ^^^«k«-k + 



2V «^ ■ 2V 

-v{k 
~2V 



-a^a-k 



w_t)(ki) +i)(ki +k2) , 
+2^ " flki+ksakiak^ 

k 

w v(ki)+u(ki + k2) ^ ^ 

+ " «ki«k2«kl+k2 

k 



1 v^' 

+ 2y ^ (5(ki + k2 - ks - k4)t)(k2 - k3)aklak2ak3ak4• 
kl,k2,k3,k4 

Here denotes the sum over all k s Note that H^{a) is the Wick 

symbol of the operator with respect to the zero mode. It is easy to compute 
also its anti-Wick symbol 

fjLf^.rrL,^. 2t)(0) 2 , ^(0) , „ W ^(0)+t)(k) 

k 

(See e.g. Appendix D for the definitions and basic properties of Wick and anti-Wick 
symbols). 

The following theorem is due to Lieb, Seiringer and Yngvason [20]. 

Theorem 3.5. Assume that the energy density e^ exists. Assume also that ■Ok is 
bounded. Then 

e„= lim y-Mnf{infspiI^(a) : a e C}. (3.13) 

Thus we can replace with the Hamiltonian H^{a) when computing the energy 
density. 

Proof. The anti-Wick symbol of the number operator with respect to the zero 
mode is 

iV^(a) = |a|2-l + ^'a^ak. 

k 

Note that for ip := supku(k), 
Now 

infspil^ < inf {infspif^(a) : a G C} 

< inf|infspF^(a) + ^iV^(a) : a e c| + ^ - ^ 

< infspiJ^ 2^ 



< infspif,^, + ^-M, 



whore in the last inequality /xi < /U and V is large enough. Dividing both sides by 
V and letting L ^ oo we obtain 

e„< lim y-Mnf{infspiJ^'(Q) : aeC}<e^,,. (3.14) 

I/— >(X) 
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Now [0,/x] 9 /xi 1-^ Cyai is a finite concave function, hence it is continuous, which 
implies (3.13). □ 



3.3. Bogoliubov method. Let us describe a version of the Bogohubov approxi- 
mation adapted to the grand-canonical approach. A similar discussion can be found 
e.g. in the review of Zagrebnov-Bru [40]. 

In what follows we will always use the grand-canonical approach. We will drop 
II from Hf^, e^(k), etc. 

For a e C, we define the displacement or Weyl operator of the zeroth mode: 

:=e-"»5+""°, (3.15) 

and the corresponding coherent vector fla := W*ri. Note that Wa is the only Weyl 
operator commuting with the momentum, and hence fla is the only coherent vector 
of momentum zero. 

Let us apply the "Bogoliubov translation" to the zero mode of H^. This means 
making the substitution 

ao = ao + cx, Uq — dQ + a, 

ak = ak, al = al, k ^ 0. (3.16) 

Note that 

and thus the operators with and without tildes satisfy the same commutation re- 
lations. In addition, the annihilation operators with tildes kill the "new vacuum" 

For notational simplicity, in what follows we drop the tildes and we obtain 

= -^|a|^ + |^|ar (3.17) 

^ , , (^(0)+^(k)) | 

[a aua-u + a a^a_^) 

Et)(k) 

k,k' 

+ 2^ '5(ki+k2-k3-k4)^^^;^ — ^«ak3ak4. 

ki,k2,k3,k4 

The expectation value of the state given by equals the constant term of (3.17), 
that is 

{n^\H^n^) = -f,\a\^ + ^\a\\ (3.18) 

(3.18) is minimized for lap = This choice kills also the linear term on the 

second line of (3.17). 
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Let us choose a so as to minimize (3.18). This means, we choose e'"^ and set 
a = e'^ . Then the Hamiltonian becomes 



2^0) + 
+ /TTTrvrr ie Ak+k'^kak' + e akak-ak+k') 

+ "^(ki + k2 - kg - k4) ' ^'^^^^ ^''^ al^al^au^au^. 

ki,k2,k3,k4 

(Note that wc have made no approximation yet). The first 3 Uncs of (3.19) form 
a quadratic Hamiltonian, which will be denoted by H^^. Now let us make the 
assumption that if^g can be treated as an approximation to H^. For a possible 
justification for this approximation see Subsection 4.1. 

It is easy to find the excitation spectrum of H^^. To this end, for k ^ we make 
the substitution 

"k = '^k^k ^ Skfo-k, Ok = Ck6k - Sk^lk) (3.20) 
Let 9 = {9u) be a sequence such that 6o = and 

^k 

Ck := cosh|6lk|, Sk := -■r^sinh|6'k|. 

Fk| 

Introduce the unitary operator 

Uo := Y[e~ ^^'''^^"■-^'^^^^'^^"■-^ . (3-21) 

k 

Note that 

U;akUe = 6k, U;alUe = bl, 
and hence ak,ak satisfy the same commutation relations as 6k, ''k- Note also that 
we have Sk = S-k and Ck = c_k = y/1 + SkSk- 

The zero mode has to be treated separately. Let us introduce the operators po 
and xo which are defined as 

= ^(e'"o^+e-'^ao), xo = -^(- e'^ < + e"'^ ao). (3.22) 

They are self-adjoint operators and satisfy the commutation relation [a;o,j?o] = i- 
As we can see they are the "momentum" and "position" of the zero'th mode. 

We choose the Bogoliubov rotation that kills double creators and annihilators, 
which amounts to 

- = ^A\^-[tJ^\ I -M ^ (3.23) 




and Ck = -\/l + |skP. We obtain 

ifbg = i^bg + MPo + 13W(k)6k^^k, (3.24) 
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where the elementary excitation spectrum is 




'knik=+2«(k)^). (3.25) 



(where the sum above includes the zero mode again). 
The lES of i7^g is given by 

ebg(k) =inf{wbg(ki) + ---+a)bg(k„) : ki + • • • + k„ = k, n=l,2,...}. (3.27) 

Note that (x'bg(k) and ebg(k) are so far restricted to ^Z'^. But these functions 
are well defined also for all values k e M"^. In fact, the thermodynamic limit of the 

lES of is simply R<^3k^ ebg(k). 
We have (in any dimension) 

(1) infk^o ^ = inf ^i(ik2 + 2^) =: Cc.,bg > 0; 

(2) limk^o "'jkl'^'' = =: Cph,bg > 0. 
Therefore, by Theorem A. 4 (1) and (2) we have 

(1) infk^^o ^^^^ = Ccr.bg; 

(2) limk^o = Cph.bg- 

Thus, ebg(k) has all the properties described in Conjecture 3.4. 
We can also compute that for small |k| 

ak«^M^/^|k|-V2. (3.28) 

(3.24) has no ground state because of the zero mode. Let ^ be any vector that 
minimizes all nonzero modes. Clearly, for k 7^ 

(*|<ak*) = |skp - ^^'^ 



2|k| 
Let 

N' := Y^ala^. 

k 

be the number of particles away from the zero mode. Clearly the density of particles 
away from the zero mode in the state ^ equals 

^mN'^) = \^Y^\s^.?. (3.29) 

k 



We expect that for large L, (3.29) converges to 



1 



(27r) 



Skfdk. (3.30) 



Note that in dimension d = 1 there is a problem with the formula (3.30), since 
|k|~-'^ is integrablc only in dimension li > 1. Therefore, for li = 1 (3.30) diverges. 
Thus, the Bogoliubov approximation is problematic for d = 1 if we keep the density 
of particles p fixed as L — > 00. To our knowledge, (3.28) and the above described 
problem of the Bogoliubov approximation in = 1 was first noticed in [8] . 

Nevertheless, in spite of the breakdown of the Bogoliubov approximation, many 
authors believe that also in d = 1 the lES exhibits the behavior e^(k) w Cph|k| with 
Cph > for low momenta, see e.g. [27], Chapter 6, [18, 17]. 
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3.4. Improving the Bogoliubov method. For '^U^ 9 k i-^ G C, a square 
summable sequence with = ^-k, let Uq be defined as in (3.21). (This time we 

allow Oq to be nonzero). For a G C, let Wa be defined as in (3.15). 

Vafi '■— UeWa is the general form of a Bogoliubov transformation commuting 
with P^. Let O denote the vacuum vector. Note that 

^a,e ■= U* Q^l 

is the general form of a squeezed vector of zero momentum. 

One of our next objectives is to look for the squeezed vector that minimizes the 
expectation value of H^. As in Section 3.3, we will also compute the Hamiltonian 
expressed in new creation and annihilation operators adapted to the new vac- 
uum ila,e- We do this in two steps. First we perform the Bogoliubov translation 
(3.16), which results in the expression (3.17). Then we perform the Bogoliubov 
rotation (3.20). This time, however, we apply it to all the modes, including 0. 

The Hamiltonian after these substitutions in the Wick ordered form equals 

= + C% + C^bo 

+ i ^ 0^(k)6£6*_k + ^ 5^ 0''(k)6k6-k + E ^''(k)6£6k 

k k k 

+ terms higher order in 6's. (3.31) 

Clearly, 

(fia,e|if^fia,e) = B^, {binc.,e\H''binc.,e) = B^ + £>^(k). 

Therefore, we obtain rigorous bounds 

E^<B^, + e^{k) < B^ + {]{.). 

If we require that B^ attains its minimum, then we will later on show that 
and 0^(k) vanish for all k. Henceforth we drop the superscript L. 



,j:(!£_,+Mi)±£MHAKP 

k ^ 

Et'(k) ,_o o_ ^ 

——(a SkCk + a SkCk) 

k '^^ 

V- w(k - k') 

+ 2V '^kSkCk'Sk' 

k,k' 

, Y^ i)(0)+t)(k-kO ,^ 2, ,2. 

+ 2^ 2v ' ' ' ' ' 

k,k' 

^ A(0)| 12 , V- (^^(0) + i^(k)) , , _ , 

C = I -^|ap -l^ + 2^ y |gk|' 1 (qco - aso) 

E^(k) / _ _ N 
— ^ (asoCkSk - acoCkSk) ; 

k 
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In order to express D{k) and 0(k), it is convenient to introduce 
/k: = y-M 

, ,^,2 t)(0)+i)(k) , i)(k'-k)+t}(0) ,^ 2 

+ y ^7 l^k'l , (-i-^^j 

k' 

nvik) v-^ wfk' — k) 
5k: = a'^-E V ^ 

k' 

(Note that /k is real). 

i?(k) = /k(4 + |5kn -CkGsk5k + SkSk), (3.34) 
0(k) ^ -2ckSk/k + 4ffk + Ck.9k. (3.35) 

The main intermediate step of the calculations leading to the above result is de- 
scribed in Appendix F. 

3.5. Conditions arising from minimization of the energy over a. We de- 
mand that B attains a minimum. To this end we first compute the derivatives with 
respect to a and a: 

OaB = |a| +2_^ |sk| ja 

v-«(k)_ 

- 2^ — SkCka, 

k 

V- m - 

- — TT^SkCkO:- 



k ^ 



Note that 



C = CodaB — SodaB, 

so that the condition 

d^B = d^B = (3.36) 

entails C = 0. The condition (3.36) yields 

^(0), ,2 sr^ viO) + v(k') . ,2 a2 '(k')_ ,^ 
H = -^|ap + ^ ^ \ ^ l ^kf - ^ E ^«k'Ck'. (3.37) 

k' k' 

This allows to eliminate fj, from the expression for /k: 



Jk : = ^ + |a' 



2 



t)(k--k)-t)(kO 2 , V-i'(k')^ ^ 

+ 2^ y l^k'l +^l.-^^Sk.Ck.. (3.38) 

3.6. Conditions arising from minimization of the energy over Sk- Comput- 
ing the derivative with respect to Sk, Sk we can use 

a „ Sk 

27 



<B = A,,-|(., + M)-,,|., (3.39) 

= A,,-|(o, + !£hJ!)-aA. (3.40) 

One can calculate that 

0(k) = f-2ck + ^) - 

Thus ds^^B = ds^B^O entails 0(k) = 0. 
(3.39) and (3.40) also imply 

Ck, 

2 

and hence 

.9kSk = Sk^k- (3-41) 
It is convenient to introduce the parameters 

^k := 2skCk, 

a := 4 + |skr. 

Now, using (3.41) we can write 

D(k) = Ck/k-5k5k, (3-42) 
0(k) = -5k/k + Ckffk. (3.43) 
Equating 0(k) to zero and assuming that /k 7^ we obtain 



Sk9s,^-B - Sk^k-B = — (ffkSk - 5kSk), 



2?(k) = sgn/k^/2 - |5kP, (3.44) 
S. - (3.45) 

a = (3.46) 

We will keep instead of fi as the parameter of the theory, hoping that one can 
later on express fi in terms of a^. We set e'"^ := j^. Then we can write 

+ ^Hi^^|l^(a,-l) + |:f).-J,„ (3.47) 

. 2v{^) ^ - k) ^ 

5k. = a-^;^-2^ — — ^-k'. (3.48) 

Then we can express fi by 

i)(0) ^;(0)+^(k') ^^ ^2^^«(k')^ ,„ 

M=— |a| +2^ (Ck'-l)-e 2^^^5k'. (3.49) 

k' k' 

One can express the minimizing conditions in the following theorem. 

Theorem 3.6. (1) Suppose that [ap > and e'"^ are fixed parameters. Let the 
first derivative of B with respect to a, a, (sk), (sk) vanish. Let /k, (?k be 
given by (3.47), (8.4.8). For any k e ^Z'' we have then f^ > Ic/kP, and 
the equations (3.44)~(3-46) hold. 
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(2) We have 



dadaB 






■ /o 


50 


dlB 


dadaB 




. 5o 


/o 



(3.50) 



(3.50) is positive/negative definite iff D{0) is positive/negative. (3.50) is 
zero iff D{0) = 0. Besides, 



^(0)=2sgn/o ma^El^^^ 



(3.51) 



In the above theorem we included all possibihtics that guarantee the stationar- 
ity of B. Clearly, the case of -D(k) < seems physically irrelevant. But this is 
equivalent to /k < 0. Therefore, under the additional condition DCk) > 0, we can 
drop sgn/k from (3.44). 

In the case of the zero momentum we have an additional argument for the posi- 
tivity of £'(0) given in (2). £'(0) > is in fact equivalent to the condition (3.50)> 0, 
which is necessary for the existence of minimum of B. 

Let us compute the ground state energy in the improved Bogoliubov method. 
Inserting (3.37) to the expression for B we obtain 



k 



2V 



k,k' 



8V 



where recall that |skP = |(Ck — !)• Using (3.37) again to eliminate |ap in favor 
of /U, and then computing the derivative with respect to fj, we obtain 

-d,B = \a\' + J2\s^\'- 



Therefore, the grand-canonical density is given by 

+ (3.52) 



3.7. Thermodynamic limit of the fixed point equation. One can ask whether 

the method described in the previous two sections has a well defined limit as i ^ oo. 
A natural way to take this limit, at least formally, involves the following steps. We 
put a = s/Vk e''^, for some fixed parameter k > having the interpretation of the 
density of the condensate. We expect Sk (and hence S'k, etc.) to converge to a 
function depending on k G M'' in a reasonable class. Finally, we replace X] by 

k 
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pip- / dk. Thus equations (3.47), (3.48) and (3.49) are replaced with 

/k = y + Kv(k) + J iv(k' - k) - «(k'))(Ck- - l)dk' 



■ 2i2nr J ^(k')^k'dk', (3.53) 

,k = /«e'2-«(k)-^^y{,(k'-k)Sk'dk', (3.54) 

^ = ^^^^^ + 2(hr J + «(k'))(Ck' - l)dk' 

/{i(k')5k'dk'. (3.55) 



2(27r 

We also obtain (in the physical case of positive D) 



D{0) = 2^v{0)n^^ I dk{,(k)5k. (3.56) 

(3.56) is typically positive - thus the quadratic part of the Hamiltonian (3.31) seems 
to have a gap. 

One can try to find a, (5k) satisfying the minimization condition by iterations. 
A natural starting point seems to be = 0. Then, by (3.37) or (3.55), /U = v{0)k. 
After one iteration we obtain 

/k = Y + KV{]S.), 

,gk = Kw(k), 
(k) = 

5k = 



£)(k) = v'(kV2)2 + k2Kt)(k), 
^^(k) 



V(kV2)2 + k2Kt)(k) 

Thus -D(k) = cjbg./j(k) given by (3.25) - we obtain the grand-canonical Bogoliubov 
approximation. 

In the case? of finite L we cannot continue iterations because of 5*0 = oo. 

In tlicrmoclynaniic limit, the value at zero may not matter, since k is a continuous 
variable, for small k behaves as ~ |k|~-'^ (this was noted already in (3.23)). In 
dimension rf = 1, if we try to do the next iteration we obtain divergent integrals. 
Thus, we cannot continue iterations. However in dimensions d > 2 the integrals 
arc convergent and we can do the next iteration (and presumably we can keep on 
going). 

The energy gap appears already at the second iteration. 

3.8. Uncorrelated states. Let Ho denote the space spanned by 1 and let W{k,-k} 
denote the space spanned by e''^^ and e"'"*^^. Clearly, 



L'{A)=no®( ® W{k,-k} • (3.57) 
V{k,-k} V 

The sum in (3.57) runs over all two-element sets of the form {k, — k} with k e ^Z''. 
The exponential property of Fock spaces yields 

r,(L2(A)) = r,(Ho)(8) f ® r,(H{k,-k})V (^-^s) 

V{k,-k} / 
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(See e.g. [35] for the definition of the tensor product of an infinite family of 
Hilbert spaces used in (3.58). Note that in each of the factors of the tensor product 
of (3.58) we distinguish a normahzed vector - the vacuum vector). 

We win say that a vector "if G rs(-L^(A)) is uncorrelated with respect to (3.58), 
or simply uncorrelated, iff it is of the form 




for some 

G r,(Ho), ^{ic-k} e r,(H{k,-k})- 

We define the uncorrelated ground state energy in the box 

E^^ := inf {(*|if^*) : is uncorrelated and of norm l} . 
For k e ^Z'' we define the uncorrelated lES in the box 

eun(k) := inf {(*|il^(k)*) - : * is uncorrelated, ||*|| = l} , 
and for k e K'' we define the uncorrelated IBS in the thermodynamic limit 

eun(k) := sup liminf inf eun(k') 

5>o\ \^k'e^Z<*, |k-k'|<5 J J 

Clearly, from the mini-max principle we obtain: 

E^<E!^^, E^+e\k)< E^^ + (k) . 

Conjecture 3.7. We believe the following statements to hold true: 

(1) The map R'' 9 k i-^ eun(k) € R zs positive and continuous away from 0. 

(2) Let k e M'^\{0}. If L ^ oo, kL G ^1'^, and ki, k, then e^„(kz,) 
eun(k). 

(3) (k) > 0. 

Thus we conjecture that using only uncorrelated states in variational determi- 
nation of excitation spectrum have serious limitations. Wc expect the results to 
be well behaved in thermodynamic limit, but they will probably not capture the 
phononic behavior at the bottom of the lES, and in particular we will obtain an 
energy gap. 

Note that the squeezed vectors fla,0 and the particle excitations over the squ- 
uezed vectors b^fla,e are examples of uncorrelated vectors. Therefore, the expecta- 
tion values of in these vectors give an upper bound on E^^ and E^^^ + e^jj(k). 
We showed that for these expectation values one should expect an energy gap - we 
expect this gap to persist even for more general uncorrelated states. 

Thus, in order to obtain more satisfactory bounds it seems that one needs to use 
corretoed vectors. Note that the Bijls-Feynman variational vector A^k^*/!! A^k\l/|| is 
correlated for k 7^ 0, even if ^ is uncorrelated. 

4. Perturbative approach 

In this section we will use the grand-canonical formalism. We replace the poten- 
tial v{x) with Xv{yi), where A is a (small) positive constant. We will drop /U from 
most symbols and instead we will make the dependence on A explicit. Thus instead 
Hjf we will write H^''^. 

31 



4.1. Perturbative approach based on the Bogoliubov method. Let us go 
back to the Bogoliubov method described in Subsection 3.3. Using the formula 
(3.19) we can split the Hamiltonian as 

2 

where 

+ Y.^^^) ^ (e"''" «ka-k + e''" ala*_^) , 

-^i •= Z^^^TKTT?'^® Ok+k'«kak'+e akflk'ttk+k')' 

ttL , 1 1 1 ^'^(k2-k3) , , 

Hi := 2^ d(ki + k2 - kg - k4) ^^^^ a^^a^^a^^^ak^. 

ki,k2,k3,k4 

Note that H^, n = — 1,0,^,1, do not depend on A. This suggests that one can 

try to apply methods of perturbation theory to compute the ground state energy 
of H^'^ treating V^Hi + XH\ as a small perturbation of the quadratic Bogoliubov 
Hamiltonian 

X-^H^^+H^. (4.2) 
It is also tempting to compute the excitation spectrum, applying perturbation 
methods to the same splitting of H^'^ restricted to the sector of fixed momentum 
k. Unfortunately, when one tries to implement this idea one encounters serious 
difficulties due to the infrared problem: the operator (4.2) does not have a ground 
state, neither globally, nor in fixed momentum sectors (because of the k = mode). 
Further on we will describe a natural approach that should help solve the infrared 
problem and should give a better starting point for the perturbation methods. 

In any case, the splitting suggests the following conjecture (which is the grand- 
canonical version of Conjecture 1.2). Let e^(k) be the grand-canonical lES for the 
potential Xv and let ebg.;j(k) be given by (3.27). 

Conjecture 4.1. Let d>2. Then for a large class of repulsive potentials we have 

JI^^mW =ebg,,.(k). 



4.2. Perturbative approach based on improved Bogoliubov method. We 

fix the size of the box, /z and e'^ and we assume that we solved the fixed point 

equation described in Subsection 3.6 and there is an energy gap. We assume that 
the solution is unique. The expression for H"^'^ Wick-ordered with respect to the 
operators 6k, allows us to write 

jjX,L ^ x-^H^i + H^'^ + \/A^^^•^ + XH^'^, (4.3) 

2 

where X~^H^'^ = B is the constant term, H^''" = X]k^(^)^k^k is the quadratic 
term, H^''" and H^''" are respectively the third and fourth order parts of H in 

2 

operators &k and sec (3.31). 

The splitting (4.3) can be used to set up a perturbatve approach for computing 
the energy density and excitation spectrum. The presence of a gap will be actually 
an advantage in this case. 
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More precisely, let us consider the following Hamiltonian 

jjX,s,L ._ _^ ^A,L _^ y/^H^'^ + SH^'^, (4.4) 

2 

where (5 is an additional parameter introduced for bookkeeping reasons. We treat 
6^^H^'^ + Hq '^ as the unperturbed operator and the rest as a perturbation de- 
pending on the small parameter 5. 5~^H^'i + Hq'^ has a ground state fia.e, and 
even a mass gap, so the perturbation expansion in terms of 5 for the ground state 
vector and energy is well defined and for small 5 

oo oo 



n=0 



where ^fg'^ = Vta,e- (It is easy to see that all powers of 6 for the energy are integral). 
At the end we will substitute A for 5: 

oo oc 

~ ^(A"*^-^ + A"+5*^;^i), E^'^ ~ X"E^'^. (4.5) 

71=0 n=—l 

Let e^(k) be the subadditive hull of £'^(k). Assume that for some ki, . . . k„ with 

k = ki H 1- k„ we have e^(k) = Z)^(ki) H \- D^{\s.n)- This implies that the 

vector (n!)~^/^6^ ' ' ' ^a.e is at the bottom of the spectrum of S~^H^'^ + H^'^ 
in the sector of momentum k. Again we can write down the perturbation expansion 
in terms of 5 for the excitation spectrum, convergent for small 5: 



M/^-^'^(k) = ^(<5"M/^-^(k) +5"+^M,^;^^^(k)), e^Ai.(k) = ^"en''^(k), 

n=0 n=—l 



where ^^'^{k) = h^Qcx,e- Then we put ^ = A obtaining 

■oo oo 

*^'^(k) ^ ^(A"*^'^(k) + A"+^*^;^ (k)), e^'^(k) ~ Yl ^"^"■''(k)- (4.6) 

n=0 n=—l 

Of course, we do not claim that the power series (4.5) and (4.6) have a nonzero 
radius of convergence. We only hope that they are in some sense asymptotic to the 
physical quantities. 

We hope that the perturbation expansions (4.5) and (4.6) survive thermody- 
namic limit. We do not expect that the nth terms of these expasions will be of 
order 0(A"). However, we hope that each next term will give a better approxima- 
tion, as expressed in the following conjecture: 

Conjecture 4.2. (1) For any n, there exist 

pX,L 

■= lim " 
e^(k) := lim e^'^(ki), k^ ^ k. 

L^oo 



(2) 



lim ell 



lime^(k) = ebg(k). 
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(3) For some < ai < (72 • • • with lim cr^ — oo, 

n — *oo 

A"e^ = 0(A''"), n=l,2,... 

A"6^(k) = 0(A'^"), n=l,2,.... 
For as above and all n, 

n 

j=o 

(1) is the existence of thermodynamic hmit of the perturbation expansion. (2) 
tells us that the lowest order terms in this expansion agree with the quantities 
obtained in the Bogoliubov approximation. (3) means that the later terms in the 
expansion are in some sense lower order than earlier. (4) says that there is no gap 
in the excitation spectrum at k = 0. 

It seems that a result similar to the above conjcctiirc could be easier to prove 
than a result about the true energy density and the true infimum of the excitation 
spectrum. 

Let us sum up the procedure that we propose to compute various quantities for 
Bose gas with A small and fixed j2. We will call it the Improved Bogoliubov Approach 

(1) Find a translation invariant squeezed state fia.e minimizing the expectation 
value of the Hamiltonian H^'^. 

(2) Split the Hamiltonian as in (4.3): 

2 

according to the power in creation/annihilation operators adapted to ^a,0- 

(3) Introduce the fictitious Hamiltonian with an additional coupling constant 
5. 

2 

(4) Compute the desired quantity perturbatively, obtaining a (formal) power 
series c^'^'^ = E„'^"c^''^- 

(5) Go to thermodynamic limit with each term in the series separately, obtain- 
ing c^^ := limL^oc c,^;-^. 

(6) Set 5 = X, obtaining the power series = X^„A"c^, which is the final 
expression for the desired quantity. 

4.3. Approach with Isolated Condensate. In the literature there are many 
works that are based on a somewhat different approach to the Bose gas with small 

A and fixed fi. This approach is sometimes called the Approach with Isolated Con- 
densate. We would like to compare it to the Improved Bogoliubov's Approach. 
Let us describe the basic steps of this approach: 

(1) Make the c-number substitution, obtaining the Hamiltonian H^'^{a) as in 
(3.12). 

(2) Substitute a = V X~^kV and split the Hamiltonian as 

H^-^{a) = X-^Hlf + H'^-^ + VXHl^ + XH^'^ . 

2 

according to the power of A. 

(3) Compute perturbatively the ground state energy, obtaining a (formal) power 
series E^^'^'^ = Y.^^ X'^E^;^. 

(4) Compute the desired quantity as a (formal) power series c^''^'^ = A"cJ^''^. 

(5) Minimize (up to the desired order in A) E^'*^'^, obtaining k^'^ as a function 
of A, L. 
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(6) Substitute k"^'^ in the expression for the desired quantity, obtaining c^'^ = 

(7) Go to thermodynamic hmit with each term of the scries separately, ob- 
taining = limi^oo c^'''^- The final expression for the desired quantity 
is 

n 

As proven by [20] (see Section 3.12), the Approach with Isolated Condensate is 
exact in thermodynamic limit for the energy density. In the case of finer quantities, 
such as the infimum of the excitation spectrum or Green's functions, we do not see 
why thermodynamic limit should make this approach exact. 

Improved Bogoliubov Approach and Approach with Isolated Condensate seem 
to have a lot in common. In both of them the main step involves calculations with a 
quadratic Hamiltonian perturbed by 3rd and 4rth order perturbations. In both ap- 
proaches the quadratic term does not contain a term linear in creation/annihilation 
operators. Note also that in both procedures the dependence of the final quantities 
on the coupling constant A can be quite complicated, and not given just by a power 
series. 

The Approach with Isolated Condensate may seem simpler technically, since the 
perturbation expansion is applied to a simpler splitting, whereas in the Improved 
Bogoliubov Approach the first step involves solving a complicated fixed point equa- 
tion. It is however quite clear, that the Improved Bogoliubov Approach is physically 
better justified than the Approach with Isolated Condensate. In the former no term 
is dropped. In the latter, at the very beginning we drop an important term from 
the Hamiltonian. 



5. Observables 

In this section we work in the grand-canonical approach. We drop the subscript 
fi and L, so that iJ^ is denoted by H. (In particular, in order not to clutter the 
notation we hide the dependence on L, which however plays an important role in 
what follows). 

5.1. Spontaneous symmetry breaking. The Hamiltonian H is invariant with 
respect to the transformation generated by the number operator e'"^^. Conse- 
quently, its ground state can be chosen to have a definite number of particles. It 
is however believed that in thcrmodynamical limit this gauge invariance is sponta- 
neously broken. (In fact, it is broken in the Bogoliubov method). 

Let us try to describe this symmetry breaking rigorously. Following Bogoliubov 
[4] , we perturb the Hamiltonian by a non-physical perturbation 

H,:=H-iyVv{a*o + ao), (5.1) 

where u > 0. depends on the gauge: 

e^^^ e-'^^ = H- uVV{e'^ a*o + e''^ ao). (5.2) 

Let us assume that has a unique ground state given by the vector ^f,^. Note 
that the Hamiltonian Hu is real, therefore we can assume to be real as well. iJ^ 
is translation invariant and the group of translations of the torus is compact. Hence 
we can take to be translation invariant. The expectation value with respect to 
the vector will be denoted 

(•). := (*.!•*.). 
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Because of the translation invariancc wc have 

(ak). - 0, k^O; 
(ak«k')i> = 0, k 7^ k'; 
(akfl-k'). = 0, k^k'. (5.3) 

Thus the nontrivial one- and two-point correlation function are 

{ao)^, = (a*)^; (5.4) 
(okOk)!/ = (a!Lka-k)i/; (5.5) 
(aka-k)i/ = (a^alk)i.. (5.6) 

and the expressions (5.4), (5.5) and (5.6) are all real. Their reality follows from the 
reality of the Hamiltonian (5.1) and the reality of 
Let us assume that there exists the limit 

(•) := hm lim (5.7) 

as a state on a suitable family 2t of observables. 

Clearly, 21 is invariant with respect to the Hermitian conjugation. Moreover, the 
group of translations e'^^ • e"''''^ and the dynamics e"'''^ • e"'''^ act on 21. 

Clearly, the ground state of e'"^^ Hi, e"^'^^ (before taking the thermodynamic 
limit) is given by e''^^^'^. Replacing "^i, with e'"^^ ^'i/ and performing the limit 
(5.7), we obtain a new state on 21. If (5.4) or (5.6) are non-zero, then this new state 
differs from : (5.4) has to be multiplied with e'"^ and (5.6) with e'^"^. 

Clearly, (5.3) are true if we replace with (•). It is natural to assume that the 
following limits exist: 

p := lim lim (5.8) 

s/li := lim lim (5.9) 
(flk^k) = lim lim {a^ak)u, k 7^ 0; (5.10) 
(akO-k) = lim lim (aka_k)j/, k 7^ 0. (S-H) 

t'XjO L—^oo 

Clearly, the expressions (5.9), (5.10) and (5.11) arc again real, p is the density and 
K can be interpreted as the density of the condensate. Both depend on p,. 



5.2. A priori estimates. We will use notation explained in an abstract setting in 

Appendix H, where the reader will also find some general remarks about Green's 
functions and their motivation. In particular, for a pair of operators A, B we the 
static Green's function is defined as 

{{A,B)), := {A{H,-E,)-'B), + {B{H,-E,)-'A),. 

Recall the operator 

A/-q := ^a*+kak = / e^''''a*axdx, 
k 



(see (3.5) and (2.12) 
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We will tacitly assume that we can perform thermodynamic limit of various 
observables, such as 



((ak,a-k)) 



lim lim {{al,ak))y, 
lim lim ((ak, a-k))i/; 



5k '= lim lim 



(^^k^^k). 



lim lim 

z/\,0 L— »oo 



(Compare with the definition of Sk and Xk in (2.15) and (2.16)). 
In this setting we have the following analog of (2.13): 



[H^, iVk]] = yiV + ^(ao + a*). 



(5.12) 



It implies the so-called f-sum rule: 



^{NZ{H^ - E,)N^), + \{N^{H, - E^)Ni)^ = y (AT), + uVV{ao).. (5.13) 
By the Schwarz inequality and taking thermodynamic limit, we obtain 



Sk < ^|k|yxir. 



(5.14) 



In the theorem below, (5.16) is due to Pitacvski and Stringari [26, 32], and (5.17) 
is the zcro-tomperature version of the famous Theorem of Bogoliubov [4] . 



Theorem 5.1. 



(OkOk) > 

> 

((ak,ak)> > 



1 



4skP 2' 

K 1 

2\k\^p ~ 2' 

K 



(5.15) 
(5.16) 
(5.17) 



+ 



((ak,a-k)) + ^ 



Proof. To simplify the presentation, our proof will be not quite rigorous, since we 
will ignore v and skip thermodynamical limit involving lim lim . 

L—>-oo 

We set A* = Ok and B := A^k in the uncertainty relation (G.2) and we obtain 



(ak«k) + :r (A^^^k)> 7l(ao)|^ 



This proves (5.15). Now (5.14) implies (5.16). 
To prove (5.17) introduce the operators 



Qk := ^k 



k5 



i?k 



i[Qk,^?]. 



We obtain 



[Qk,[H,Qk]] = 2k'N-k'Q2k, 
[<3k,ak] = — flo — a2k, 
[Qk,a*_u] = ao + al2k- 
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Therefore, 



^((i?k,i?k)) = ^([Qk,[i?,gk]) 



= (iV)k2, (5.18) 

((ak,-Rk)) = i([Qk,ak]) 

= -i(ao>, (5.19) 

((a*_k,iik)) = i([Qk,a*_k]) 

= 1(4) • (5-20) 

(5.19) and (5.20) arc sometimes called Bogoliubov sum rules [4, 32]. 
For a complex parameter t we have 

-i(ao)+it(aS) = (((ak + tolk).^k)). (5.21) 
We take the square of the absolute value of (5.21), apply (G.3), and we obtain 

\{a^)\'~t{a^r-t{aor + \fr\{ao)\' 

< (((flk + talk), « + ta-k)))((i?k, Ru))- 
Taking into account (5.18), wo obtain 

l(«o)P 



< ((ak,a£)) 



+M ((«-k,ak)> + 



(Ar)k2 

(A^)k2 ; 



+*((K,a-k)) + ^) 



Using 



(ao) = {a*Q), ((a^k, o_k)) = {{al, Ok)), ((a^k, O) = (K, a-k)), 
we obtain 



((ak,a_k)) + 



(ao)' 



(Af)k2 



< ((ak,ak)) - 



l(ao)| 



2 



(iV)k2 ' 

which implies (5.17). □ 
Note the following consequence of (5.17): 

((a^,ak))-((a^,alk)> > (5.22) 

Theorem 5.2. Let e(k) 6e i/ie /iSS' at momentum k. T/ien 

e(k)' < ^(y-'^ + MO) (5.23) 

"^2(^0^ / ^^^^ (2(<+kaq+k) + (<+ka-q-k) + (oq+ka-q-k)) dq), 

e(k)2 < (^^^)' + 2k2yM!(«;«^)dq (5.24) 
+p J dx(l — coskx)V?^t'(x)(aQa*axao), 



direction o/k. 



where k denotes |k| ^k and V?'w(x) denotes the second derivative of v in the 
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Proof. To prove (5.23), we use (G.8) with A* = Ok — alk and B = A^k and then go 
to thermodynamic Kmit. 

To prove (5.24) wc use (G.9) with A = N^^ □ 
Both estimates of Theorem (5.2) indicate the phononic character of the excitation 
spectrum. The estimate (5.23) is due to Wagner [37, 32]) and involves the symmetry 
breaking parameter k. The estimate (5.24) involves the kinetic energy and the pair 
correlation function (a^a'^a-x^ao) (here refers to the position), but does not involve 
K, hence it can be applied to situations without symmetry breaking. This estimate 
comes from [28, 24], sec also [32]. 

5.3. Green's functions of the Bose gas. Let us consider two-point Green's 
functions of the Bose gas. We assume that (•) is the state obtained by the limiting 
procedure in thermodynamic limit, and we will ignore the complications due to 
thermodynamic limit. 

We define a 2 x 2 matrix of Green's functions 

" Gii(^,k) G2i(^,k) " 

Gl2{z,k) G22(2,k) 



G{z,k) :-- 



G2i(2,k) 
Gi2(^,k) 

G22(z,k) 



= {akiH -E- z)-\l) + {al{H -E + z)-^ak), 

= {a\{H - E - z)-'al) + {aUH -E + z)-'a\), 

= (ak(i? -E- z)-ia_k) + {«-k(i? -E + z)-^a^), 

= {a*_^{H -E- z)-ia_k) + (a_k(ir -E + z)-^aU 



Note that, using the notation of Appendix H, 

G,,(z,k) = Ga„b,(^), 



where Ax := Ok, A2 



and Bi :— at, B2 :— a_ 



We use the conventions 



described in this appendix for the meaning of Green's functions both away from 
the real line and on the real line. 

It is a general fact, which does not depend on the details of the system, that 



Gii(^,k) =Gn(^,k) =G22(-^,-k), 

Gi2(^,k) =G2i(0,k) =Gi2(-^,-k). 

By the reflection invariance of the Bose gas 

Gi, (^,k) = Gij{z, -]<.). (5.25) 

Obviously, for any observable A, (A*) = (A). But the state (•) is real, hence 
(A) = (A). Note also that H — H , a\s_ = a\s_. Therefore, 

Gi2(z,k) = G2l(z,k). 

Note that Gii(0,k) = ((a^,ak)) and Gi2(0,k) = ((ak,a-k)), hence by (5.22) 



Gii(0,k)-Gi2(0,k)> 



k2- 



(5.26) 



Let us introduce the "full mass operator" 
S(z,k) 



5]ii(z,k) Ei2(z,k) 
E2i(^,k) S22(2,k) 



1 

2^ 



= (Gii(^,k)G22(^,k) - Gi2(^,k)G2i(^,k))-' 
Consequently, 

5]ii(0,k)-Ei2(0,k) = 



Gn(z,k) Gi2(z,k) 

G2l(^,k) G22(-2,k) 

G22(2,k) 



-G2l(^,k) 



-Gi2(^,k) 

Gn(^,k) 



27r(Gii(0,k)-Gi2(0,k))- 
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(5.26) implies 

— >Sii(0,k)-Si2(0,k)>0, 

ZTTC 

and in particular 

Eii(0,0)-Ei2(0,0) = 0. (5.27) 

(5.27) was first proven in the framework of perturbation theory for the Bose gas 
with isolated condensate, and is sometimes called the Hugenholz-Pines Theorem, 
[13], sec also [8]. The proof that we present is valid for the correct Hamiltonian of 
the Bose gas and is due to Bogoliubov [4]. 

Note that (5.27) implies that G{z,\c) has a singularity at (^;,k) = (0,0), which 
is an argument for the absence of a gap in the excitation spectrum. Bogoliubov 
[4] gives also an argument for the phononic shape of the excitation spectrum. The 
argument is based on the assumption that S(^,k) is regular in ^, k around (0,0). 
Note that 

detl](2:,k) = I]n(^,k)E22(^,k)-Ei2(z,k)S2i(^,k) 

is invariant with respect to the transformations k — k and z i— > —z. Finally, by 
(5.27), we know that S(0, 0) = 0. Therefore, 

detS(.2,k) = +/3k2 +0(1^1^ + |k|^). 

We have detS(^;, k) = detS(;2, k). Hence 7 and /3 are real as well. For purely 
imaginary nonzero z, det'S{z,\s.) is nonzero. Hence 7 and /3 cannot have the same 
sign. Therefore, S = > 0. 

Assume now that (3, 7 are not zero. Then < 6 < 00, and 

detS(\/^|k|,k) = 0(|k|^). 

Hence, for small w, k, the Green's function G{u), k) has a sharp peak along a'(k) = 

Appendix A. Energy-momentum spectrum of quadratic Hamiltonians 
Suppose that we consider a quantum system described by the Hamiltonian 

H = [ a;(k)a;;;akdk, (A.l) 

with the the total momentum 

P = kakflkdk, 

both acting on the Fock space rg(L^(R'')). We will call the function ui appearing in 
H the elementary excitation spectrum of our quantum system and we will assume 
it to be nonncgativc. 

Clearly, the ground state energy of H is 0. The excitation spectrum of (A.l) 
is not arbitrary - it has to be a subadditive function. This appendix describes a 
number of easy results about subadditive functions. They arc quite straightforward 
and probably they mostly belong to the folk wisdom. However, we have never seen 
them explicitly described in the literature, and we believe them to be relevant for 
physical properties of Bose gas. 

The Hamiltonian of interacting Bose gas is not purely quadratic. Nevertheless, 
some arguments indicate that the infimum of its excitation spectrum in thermo- 
dynamic limit is subadditive. A heuristic argument in favor of this conjecture is 
described in the next appendix. 

Even if one questions this argument, there exists another motivation for a study 
of subadditive functions. Quadratic Hamiltonians are often used in statistical 
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physics as approximate effective Hamiltonians. In particular, this is the case of 
the Bogohubov Hamiltonian Hhg. 

We will show that there exists a large class of subadditive functions with the 
properties properties described by Conjecture 1.1 or 3.4, (3) or (3)', and (4) (which 
correspond to the superfluidity or periodicity, and to a finite speed of sound). We 
will also show that if the elementary excitations possess these properties, then so 
does the lES. 

Let K'' 9 k I— > e(k) e M be a nonnegative function. We say that it is subadditive 

iff 

e(ki+k2)<e(ki) + e(k2), ki,k2GM'^. 

Let R** 9 k I— > a;(k) s R be another nonnegative function. We define the 
subbadditive hull ofu to be 

e(k) := inf{w(ki) + • • • + a;(k„) : ki + • • • + k„ = k, n = 1, 2, . . . }. 

Clearly, e(k) is subadditive and satisfies e(k) < uj(k). 

Clearly, if a;(k) is the elementary excitation spectrum of a quadratic Hamiltonian, 
and e(k) its subadditive hull, then e(k) is the infimum of its excitation spectrum. 

Let us state and prove some facts about subadditive functions and subadditive 
hulls, which seem to be relevant for the homogeneous Bosc gas. 

Theorem A.l. Let f be an increasing concave function on [0,oo[ with /(O) > 0. 
Then /(|k|) is subadditive. 



Proof. 



/(|kl+k2|) < /(|ki| + |k2|) 



< /(|ki|) + /(|k2|). 

□ 

We can generalize Theorem A.l to periodic functions. 

Theorem A. 2. Let f be an increasing concave function on [0, ^] with /(O) > 0. 
Define e to be the function on periodic with respect to the lattice such that if 
k e [— |, U"^, then e(k) = /(|k|) (which defines e uniquely). Then e is subadditive. 

Proof. We can extend / to a concave increasing function defined on [0, oo[, e.g. by 
putting fit) = f{% for t > ^. 

Let ki,k2 e W^. Let pi,p2 e [-|, \Y such that k^ - G Z''. Let p G [-\, \Y 
such that kj + k2 - p e U^. Note that |p| < |pi + P2|. Now 

e(ki+k2) = /(|p|)</(|pi+P2|) 

< • • • 

< /(|pi|) + /(iP2|) = e(ki) + e(k2), 

where in . . . we repeat the estimate of the proof of Theorem A.l. □ 

Obviously, we have 

Theorem A. 3. Let eo be subadditive and < uj. Let e be the subadditive hull of 
CO. Then eo < e. 

In the case of the Rose gas with repulsive interactions we expect that the ex- 
citation spectrum may have resemble that of a quadratic Hamiltonian with the 
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properties described by the following two theorems, which easily follow from The- 
orems A.l, A.2 and A.3: 



Theorem A. 4. Suppose that w >0 is a spherically symmetric function on M** and 

e is its subadditive hull. 

(1) e is spherically symmetric. 

(2) If inf ^4tt^ = c, then inf ^itt^ = c. 

■' k^^iO k/o 1*^1 

(3) //liminfk^o = c, then e(k) < c|k|. 

(4) Suppose that for some c > 0, we have ti;(k) > cmin (|k|, 1). Then 

liminf } = Cr,h implies lim = Coh- 
k^O |k| k^O |k| 

Theorem A. 5. Suppose that iv >0 is an even function on M. periodic with respect 

to Z. Let e he its subadditive hull. 

(1) e(k) is even and periodic with respect to Z. 

(2) If, for some c > 0, we have a;(k) > cdist(k, Z), then 

liminf —^—^ = Cr,h implies lim ^ = Cr,h- 
k^O |k| k^O |k| 

Appendix B. Subadditivity of the excitation spectrum of interacting 

BOSE GAS 

In this appendix we describe a heuristic argument in favor of Conjecture 3.4 (5). 
Recall that this conjecture says that the lES of interacting Bose gas in thermody- 
namic limit should be subadditive. Clearly, this would be true if the Bose gas was 
described by a quadratic Hamiltonian of a form A.l. Wc will see, however, that 
this conjecture follows as well from an assumption saying that one can describe 
excitations by approximately localized operators. 

Consider Bose gas in a box of side length L where L is very large. Let $0 
be the ground state of the Hamiltonian and Eq its ground state energy, so that 
i?$o = Eq^o and P$o = 0. Let (£^0 + ei,ki) e sp(iJ,P), i = 1,2. We can find 
eigenvectors with these eigenvalues, that is, vectors $i satisfying H^i = (Eo+ei)^i, 
P^i = ki<i>i. Let us make the assumption that it is possible to find operators Ai, 
which are polynomials in creation and annihilation operator smeared with functions 
well localized in configuration space such that PAi « Ai{P + k,), and which ap- 
proximately create the vectors from the ground state, that is w Ai<pQ. (Note 
that here a large size of L plays a role). By replacing <J>2 with e'^^ $2 for some y 
and A2 with e'^^ A2 e^'^^, we can make sure that the regions of localization of Ai 
and A2 are separated by a large distance. 

Now consider the vector $12 := AiA2^o- Clearly, 

P$12 « (ki +k2)$12. 

$12 looks like the vector $j in the region of localization of Aj, elsewhere it looks 
like $0 • The Hamiltonian H involves only expressions of short range (the potential 
decays in space). Therefore, we expect that 

H^i2 « (£'0 + ei + e2)$i2. 

If this were the case, it would imply that (£^0 +61+ 62, ki + k2) G sp(iJ, P), and 
hence would show that the lES is subadditive. 

Clearly, the argument we presented has its weak points - it is based on approxi- 
mate locality, which can be violated because of correlations due to the Bose-Einstein 
condensation. Nevertheless, we have the impression that many physicists believe 
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that even in the interacting case, in the thermodynamic hmit, one can often "com- 
pose excitations" in a sense similar to the one described above. (See the discussion 
of the concept of elementary excitations in interacting systems by Lieb [17]). 



Appendix C. Speed of sound at zero temperature 
It is well known (e.g. [11]) that at any temperature the speed of sound is given 




by 



where p is the pressure, p is the density and S is the entropy. (Recall that we 
assume that the mass of an individual particle is 1). 

Let E denote the ground state energy (which corresponds to the total energy at 
zero temperature), V the volume and n the number of particles. Note that n = Vp 
and E ~ Ve{p), where e(p) denotes the energy density. At zero temperature the 
pressure is given by 

dE I 



Clearly, at zero temperature the entropy is zero. Therefore, 

dp 

dp s=o dp 



T=0 



= A(_e(p)+pe'(p))=pe"(p). 

Appendix D. Wick and anti-Wick symbol 

Let a*, 02, . . . , a* and ai, a2, . . . ,an be creation/annihilation operators. Let H 
be an operator given as a polynomial in these operators. We can write H in two 
ways: 

(We use here the multiindex notation, e.g. a" = a"^ ■ ■ ■ a^"). Then the function 

C" 9 a H{a) = h^^soPa^ 

7,(5 

is called the Wick symbol of the operator H. (Synonyms: lower symbol, normal 
symbol, covariant symbol, a*,a-symbol, Q-representation) . The function 

7,(5 

is called the anti-Wick symbol of the operator H . (Synonyms: upper symbol, anti- 
normal symbol, contravariant symbol, a,a*-symbol, P-represcntation). 
Introduce the standard coherent states: 

W{a) := exp [ ^(— a^a* + aiOi) j , f2a := W{a)Q.. 

Note the identities (that can be used as alternative definitions of the Wick and 
anti-Wick symbols): 

H{a) = {n^\Hn^), (D.l) 
H = [ H{a)\n„){n„\ — . (D.2) 
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Let H he a, bounded from below self-adjoint operator. We have the following 
lower and upper bound for the ground state energy of H, which follow immediately 
from (D.l) and (D.2): 



mi{H{a) : a e C"} < inf spil < inf{ff(a) : a e C"}. 

Appendix E. Bogoliubov transformations 

In this appendix we recall the well-known properties of Bogoliubov transforma- 
tions and squeezed vectors. For simplicity we restrict ourselves to one degree of 
freedom. 

Let a*, a are creation and annihilation operators and O the vacuum vector. 

Recall that [a, a*] = 1 and ail = 0. 

Here are the basic identities for Bogoliubov translations and coherent vectors. 
Let 

W„ := e-«''*+«" . 

Then 

WaaW* =a + a, 
Waa*W* =a*+a, 

W*n =e-^e«''*0. 
Here are the basic identities for Bogoliubov rotations and squeezed vectors. Let 

Ue := e-f«*°*+f««. 

Then 

UeaUg = cosh |6'|a+ ilysinh 16*10*, 
U0a*U^ = cosh|6'|a* + ||| sinh|6'|a, 

U^n = (1 + tanh^ \e\)i e"^ 1^1"*"* O. 

Vectors obtained by acting with both Bogoliubov translation and rotation will 
be also called squeezed vectors. 

Appendix F. Computations of the Bogoliubov rotation 

In this appendix we give the computations of the rotated terms in the Hamil- 
tonian used in Section 3.4. 



Uoal^akUg = |sk 
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+cialak - CkSka^alk - CkSkOkfl-k + \sk\'^ a*_f^a-k, 
Ueal-at^Ug = -SkCk 

+c^a^a*_y^ — CkSkflkflk — CkSkfxlk'^-k + s^a-kCik; 
Uea\^a-\JJQ = -SkCk 

-t-Ckakfl-k — CkSkflkak — CkSkolkC'-k + s\a*_]^a^; 



t^eak+k'«kak'C^| = (co(|skP'5(k') + |sk' p(5(k)) + soCkSk(5(k + k')) oq 

- (so(|sk|'(5(k') + |skf <5(k)) + coCkSk<5(k + k')) a* 

+higher order terms; 
C/9ak+k'a£a£'?70 = (co(|sk|^<5(k') + |sk'|^(5(k)) + soCkSk<5(k + k')) 

~ (so(|skr(5(k') + |skf 5(k)) + coCkSk<5(k + k')) ao 

+higher order terms; 

(5(ki + k2 - ks - k4)Z7eakiak2"k3ak4C^0 
= CkiSkiCk3Sk3(5(ki + k2)(5(k3 + k4) 

+ kkj'iSk.r (<5(kl - k3)5(k2 - k4) +(5(ki - k4)5(k2 - ks)) 

+ (^SkiCki(sk3Ck3alk3a-k3 - Sk3«k3a-k3 - Ck3ak3a-k3 + Sk^Ck^al^ak^) 

+Sk3Ck3(skiCkiaki«Jki -Ski^kifl-ki - Ck^ak^^a^lk^ + Ski Cki a^lk^ fl-ki 
x<5(ki +k2)5(k3 +k4) 

+ (iSkal^dCkil^dki^ki - CkiSkiOkia^lki " CkiSkiakiO-ki + |Ski|^alki<^-ki) 
+ |Ski|^(|Ck2|^ak2"k2 - CkjSksCtka'^-ks " CkjSkjaksa-ka + |Sk2pa-k2«-k2)) 
X (5(ki - k3)^(k2 - k4) + <5(ki - k4)^(k2 - ks)) 

+higher order terms. 

Appendix G. Operator inequalities 

Let us fix a vector * and let (A) denote Let [A,B]+ := AB + BA 

denote the anticommutator. Occasionally, we will write [A, := AB — BA for 
the usual commutator. 

Theorem G.l. Suppose that A, B are operators. We have the following inequali- 
ties: 

(1) Schwarz inequality for an anticommutator 

\{[A%B]+)\' <{[A\A]+){[B\B]+). (G.l) 

(2) Uncertainty relation for a pair of operators 

\{[A\B])\''<{[A,A*]+){[B,B*]+). (G.2) 
Proof. We add the inequalities 

< (A + tB)*{A + tB), < ± tB){A ± tB)*, 

obtaining 

< [A,A*]++i\B\A]±+t[A\B]± + \t\'^[B\B] + . 

Then we take the expectation value of both sides and set t = — • D 

Suppose that H is an operator bounded from below and ^ is its ground state 
vector: 

= E^, H-E>0. 
Assume that (A) = (B) — 0. Then we will write 

{{A, B)) := {A{H - E)-^B) + {B{H - E)-^A). 

45 



Theorem G.2. 

\{{A*,B))\' < {{A*,A)){{B*,B)); (G.3) 
\{[A\B])\^ < {{A*,A)){[B\[H,B]]). (G.4) 

Proof. To sec (G.3) wc apply the Schwarz inequality to the positive definite form 

{{A*',B)). 

To obtain (G.4) we first use the identity 

{[A*,B]) = {{A*,[H,B])), (G.5) 

and then (G.3). □ 

Theorem G.3. Let Ho be the space 

Ho := {h{H)A^ + g{H)A*^ : f,gr' 

(the smallest invariant subspace of the operator H containing A'^ and A*^). Let 
e :=infspi? - E. Then 

T~Co 

{[A*,[H,A]]) 

< <[f;;[^;f]]), (G.7) 

{{A*, A)) ' ^ ' 

2 ([A*.[//..l]])([ir.[J/.L-]]) 

' - — kr%wp — ' ^^-'^ 

, {[[A\H][H, [H,A]]]) 

' - ^\A%{B.A\\) ■ (^-^^ 
Proof. To prove (G.6) we add 

e{A*A) < {A*{H - E)A), e{AA*) < {A{H - E)A*). 
To prove (G.7) wc add 

e^{A*{H - E)-^A) < {A*{H - E)A), e^{A{H - E)-^A*) < {A{H - E)A*). 

(G.8) follows from (G.7) and (G.4). □ 
(G.6) is called the Feynman bound [32] and (G.8) is due to Wagner [37, 32]. 
(G.9) comes from [28, 24]. 

Appendix H. Green's functions 

We consider a quantum system described by a bounded from below Hamiltonian 
H. We assume that it has a unique ground state (*| • ^f) = (•), = E'^. If A is 
an operator, then we will write 

A{t) :=e"^Ae-"^. 

The time- dependent Green's function associated to a pair of operators A, B is 
defined as the function depending on t e M 

Gtsit) = e{-t){A{0)B{t))+e{t){B{t)A{0)) 

= 9i-t){Ae'''-"-'^^ B) +e{t){Be-'*^"-^^ A), (H.l) 

where is the Heaviside function. 

We also introduce the energy-dependent Green's function, which is the Fourier 
transform of (H.l) (with one of conventional normalizations). It is the distribution 

46 



on a; e M defined as 



Ga,b{^) = limiy"GA,BWe-"^*-^l*ldi (H.2) 

= lim i y" (^{A Q-itiH-E-u^-ie) ^ ^5 ^it{H-E+u;-ie) 

= (A{H - E-uj- iOy^B) + {B{H - E + uj - iO)"^^). 

Finally, the analytic Green's function is defined for 2; G C\ (sp(-ff — i?) U sp(iJ — H)) 
and is defined as 

GaA^) = {A{H-E-z)-^B) + {B{H-E + z)-^A). 

The distribution GA,B{to) is the boundary value of the analytic function G'^g{z), 
provided that we approach the the real line from the appropriate side. Besides, in 
the energy gap both functions coincide: 

GA,Bi(^) = GXb(^-^Q), co€sp{E- H) (z]-oo,0]; 
Ga,bH - G^bH, ujeR\{sp{E-H)Usp{E-H)}; 
Ga,b{(^) = GXB(t-' + iO), Lo esp{H - E) c]0, 00]. 

Motivated by the above relations, following the usual convention, we can treat 
Ga,b{'^) and G'^g{z) as restrictions of a single fuction and drop the subscript an. 
Note that 



Ga,b{z) = Gb*,A'{z) = Gb,a{-z). 

Green's functions are well motivated physically. Let us briefly describe their two 
separate physical applications. 

Following [4], let us first describe the physical meaning of the static Green's 
function 

Gab(0) = ((AS». 

Suppose that ^ is an eigenvector of H (not necessarily a ground state). Let B 
be a perturbation with 1{e}{H)B'^ = 0. Suppose that it is possible to apply 
perturbation theory to the family := H + tB obtaining an analytic familly 
of eigenvectors \E'x with eigenvalues E^- such that Eq = E and \l/o = The 
Rayleigh-Schrodinger perturbation theory says that 

^^ = ^f + r(iJ-£;)-^B* + 0(T2). (H.3) 

Let {A)t := {'^t\A'^t)- (H.3) implies that for any operator A we have 

Thus {{A, B)) measures the linear response of eigenvalues of a quantum system. 

Let us describe a typical illustration of the physical meaning of Ga,b{z) for a 
general z. Suppose that at time the system described by a Hamiltonian H is in its 
ground state. We perturb the Hamiltonian by a weak perturbation \B and at time 
t we measure the observable A. The shift of the expectation of the measurement is 

5^{t) := (e'*(«+^^)^e-«(^+^^))-(^) 

« iA f M[B{u),A{t)]) 
Jo 

= iA y * ds ((B e-(^-^) A) - {A e--(^-^) B}) , 
47 



where we took the leading term in A. For some Imz < 0, we compute the Laplace 
transform of Sx{t), make the linear approximation and change the variable u = t—s: 



/•OO 



e-'''Sxit)dt « iAy"°°e-'*"d<y'((Se'"(^-^)A)-(Ae-'"(^-^)B))ds 

/■CX3 /-OO 

= iA y e-™^ du J (^{B e''^"-^-'^ A) - {^Ae-''^"-^+'^ B)j ds 
z 

Thus, Ga,b{z) measures the linear response of the dynamics a quantum system. 

H.l. The van Hove formfactor. A typical experiment measuring excitation 

spectrum involves scattering with a beam of particles, typically neutrons. Following 
van Hove [36] let us try to describe such an experiment mathematically. 

We can assume that a neutron of mass m interacts with each particle of the Bose 
gas through a potential w. Its incident momentum is pi. We measure scattered 
neutrons of momentum pf 7^ pi. The space describing the Bose gas and a single 
neutron is rs(A) (g) L^{M.'^) and the Hamiltonian is 

= 1+1(8) + \ j a;axw(x - y)dx, 

where y denotes the position of the neutron, Dy its momentum, A is small, and as 
usual H is the Hamiltonian of the Bose gas. Let i>p denote the plane wave function 
of momentum p, that is $p(y) = e'P''. 

Suppose that the initial state of the composite system is ^' (g , where ^ is the 
ground state. Let E be the ground state energy and Ui = the energy of the 

incident neutron. After time 2T the evolved state is given by 

e(T,pi) = e-'^"^^^ * O 

w e"'^"^^" * O $pi 

j-2T 

-iA / e-'-^TH,+-n(H,-E-a,) J^^^ ^ 
Jo 

= e-^2^(^+'^') * ® $p, 

Hq — E — Ui 

where we used the so-called Born approximation. Let Uf := j^Pf be the final 
energy of the neutron. We introduce also the momentum and energy transfer 

q = Pi - Pf , w := CTi - (Tf . 

To obtain the amplitude of the measurement of the momentum pf we take the par- 
tial scalar product of Q(T, pi) e rs(A) O L'^{R'^) with with $pj e obtaining 
the vector in rs(L^(A)) equal 

e(T,a;,q) := ($p, |e(T, p;)) 

= _ 2Ai ^_ir(H„+ir+.,+.o sinr(go -E-u:) 1 f _ 

V Ho-E-uj J J ^ 
^ _2Ai ^_,^(^,^^^,,^,,) sinT(go - E- ^) ^^^^^^^_ 

V Hq — E — uj 
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Note that the number of states in a cube dqi • • • dq^ equals equals V{2n) ''dqi • • • dq^. 
Therefore, the scattering crosssection per unit time in the Born approximation is 

^||e(r,c^,q)||V(2^)-'^ 



(N) 



1/2 

(27r)i-'^AV|^^'(q)|'5(a;,q), (H.5) 



where p = -^yJ- is as usiial the density and 

5'(w,q) = (iV)-i(*|iV*5(i?-£;-w)iVq*) (H.6) 

is sometimes called the van Hove formfactor. 

It is interesting to note that (H.5) depends on the incoming and outgoing data 
only through the momentum and energy transfer. 
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